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Interactions induced by electromagnetic fluctuations, such as van der Waals and Casimir 
forces, are of universal nature present at any length scale between any types of systems 
with finite dimensions. Such interactions are important not only for the fundamental 
science of materials behavior, but also for the design and improvement of micro- and 
nano-structured devices. In the past decade, many new materials have become available, 
which has stimulated the need of understanding their dispersive interactions. The field 
of van der Waals and Casimir forces has experienced an impetus in terms of developing 
novel theoretical and computational methods to provide new insights in related phe¬ 
nomena. The understanding of such forces has far reaching consequences as it bridges 
concepts in materials, atomic and molecular physics, condensed matter physics, high en¬ 
ergy physics, chemistry and biology. In this review, we summarize major breakthroughs 
and emphasize the common origin of van der Waals and Casimir interactions. We ex¬ 
amine progress related to novel ah initio modeling approaches and their application in 
various systems, interactions in materials with Dirac-like spectra, force manipulations 
through nontrivial boundary conditions, and applications of van der Waals forces in 
organic and biological matter. The outlook of the review is to give the scientific com¬ 
munity a materials perspective of van der Waals and Casimir phenomena and stimulate 
the development of experimental techniques and applications. 


PACS numbers: 34.20.Cf,78.67.Uh,81.05.ue 
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I. INTRODUCTION 


Phenomena originating from electromagnetic fluctuations play an important role in many parts of science and 


technology. The Casimir effect, first predicted as an attractive force between neutral perfect metals (Casimir 19481, 


has made an especially large impact. This non-classical electromagnetic force is typically associated with the coupling 
between objects with macroscopic dimensions. The same type of interaction known as a Casimir-Polder force concerns 


atom/surface configurations (Casimir and Polder 1948). The conceptual realization of the Casimir and Casimir-Polder 


effects, however, is much more general. The connection of such interactions with broader definitions of ’’dispersion 
forces” establishes a close relationship with the van der Waals (vdW) force ( Barton[ 1999[[Mahanty and Ninhain] 1976 


Parsegian, 2006). The common origin of vdW and Casimir interactions is directly related to their fluctuations nature. 


since at thermodynamic equilibrium the electromagnetic energy of dipoles (associated with the vdW force) can also 
be associated with the energy stored in the electromagnetic fields (the Casimir regime), as illustrated schematically 
in Fig. [T] The point that these constitute the same phenomenon was realized by several authors, including (Barash 
and Ginsburg 1984) who write, ’’The fluctuation nature of the van der Waals forces for macroscopic objects is largely 


the same as for individual atoms and molecules. The macroscopic and microscopic aspects of the theory of the van 
der Waals forces are therefore intimately related.” 

This ubiquitous force, present between any types of objects, has tremendous consequences in our understanding 
of interactions and stability of materials of different kinds, as well as in the operation of devices at the micro- and 
nano-scales. The Casimir force becomes appreciable for experimental detection at sub-micron separations. This is 
especially relevant for nano- and micro-mechanical devices, such as most electronic gadgets we use everyday, where 


stiction and adhesion appear as parasitic effects ( 

Buks and Roukes 

to 

o 

o 

)—‘ 

). The Gasimir force, on the other hand, can 

be used to actuate components of small devices without contact (( 

”han et al. 

2001 

)• 


vdW interactions are recognized to play a dominant role in the stability and functionality of materials with chem¬ 
ically inert components, especially at reduced dimensions. The most interesting recent example has been graphene 


and its related nanostructures (Novoselov et al. 2004). The graphene Dirac-like spectrum together with the reduced 


dimensionality are responsible for novel behaviors in their Casimir/vdW forces. The graphene explosion in science 
and technology has stimulated discoveries of other surface materials, including 2D dichalcogenides, 2D oxides or other 
honeycomb layers, such as silicene, germanene, or stanene, where dispersive forces are of primary importance. En¬ 
gineering heterostructures with stacking different types of layers is an emerging field with technological applications 


via vdW assembly (Geim and Grigorieva 2013). Other materials with Dirac spectra are also being investigated. For 


example, topological insulators, Ghern insulators, and Weyl semimetals are very interesting for the Gasimir/vdW field 
as the surface of such materials has a distinct nature from the bulk. 

The importance of the vdW interaction extends to organic and biological matter. Perhaps the adhesion of the 
Gecko, a type of lizzard from the Gekkota infraorder, has become a pop-cultural poster child for such interactions 


after Autumn and coworkers (Autumn et al. 


2002) in a series of experiments showed that complex hierarchical nano¬ 


morphology of the gecko’s toe pads (Lee 2014) allows them to adhere to hydrophobic substrates (Autumn and Peattie 


2002). Dispersion forces play an important role in the organization of other bio-systems, such as cellulose, lignin. 


and proteins. The stability of biological matter via an array of lipid membranes coupled through the vdW force is a 
fundamental problem of much current interest in soft matter physics. 

The stability of many hard materials, including composites and heterostructures, is also closely related to their 
Gasimir and vdW interactions. The electromagnetic nature makes this phenomenon inherently long-ranged as it 



FIG. 1 (Color Online) Schematic representation of dispersive interactions induced by electromagnetic fluctuations for: (a) the 
dipolar vdW force between atoms and molecules; (b) the Casimir-Polder force between atoms and large objects; and (c) the 
Casimir force between large objects. For small enough separations one can neglect retardation effects due to the finite speed 
of light c, which corresponds to the vdW regime. For large enough separations, retardation effects become important, which is 
charactersitic for the Casimir regime. 
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depends in a complicated manner upon the electromagnetic boundary conditions and response properties of the 
materials. Metallic and dielectric structures of nontrivial shapes lend themselves as a platform where this aspect 
can be investigated in order to tailor this force in terms of its magnitude and sign. The payoff is highly beneficial 
in the context of being able to reduce the unwanted stiction and adhesion in nano electro-mechanical and micro 
electro-mechanical devices and improve their performance. Structured materials, including metamaterials, photonic 
crystals, and plasmonic nanostructures, on the other hand, allow the engineering of the optical density of states and 
magnetic response, which is also useful for manipulating the Casimir force. 

Research published in the past decade has shown that the role of materials can hardly be overestimated when it 
comes to the description and understanding of dispersive interactions. In addition to recent books discussing basic 


concepts (Bordag et al. 2009b Buhmann, 2012a|b Dalvit et aL\ 2011a Parsegian 2006 Simpson et al.. 2015), there 
are several existing topical reviews on the Casimir effect with different emphasis. Aspects such as the quantum 


and Welsch 2007), experimental progress (Lamoreaux 


(PFA) with related experiments (Klimchitskaya et al. 


field theory nature (Bordag et al. 2001 Milton 2004), the quantum electrodynamics (QED) method (Buhmann 


Dalvit et al. 2011b Reid et al. 2013a Rodriguez et al. 2011a, 2014) have been summarized. 


20051, the Lifshitz theory and proximity force approximation 


2009), and non-trivial boundary conditions (Buhmann 2012b 


Nevertheless, the materials perspective of vdW/Casimir interactions has not been considered so far. With recent 
advances in materials science, especially in novel low-dimensional materials, composites, and biosystems, this held 
has become a platform for bridging not only distance scales, but also concepts from condensed matter, high energy, 
and computational physics. There is an apparent need for discussing progress beyond the existing topical reviews via 
a materials perspective and give a broader visibility of this held. The purpose of this article is to summarize advances 
in the development and application of theoretical and computational techniques for the description of Casimir and 
vdW interactions guided and motivated by progress in materials discoveries. Each section of this review describes a 
separate direction dehned by the type of systems, length scales, and applications of vdW/Casimir phenomena. An 
integral part is a succinct presentation of hrst principles and coarse grained computational methods highlighting how 
the distance scale is interconnected with adequate micro and macroscopic description of the materials themselves. 
Our intention is to stay within the equilibrium conditions and not include thermal non-equilibrium vdW/Casimir 
effects, nor critical Casimir interactions. The complexity of these omitted aspects of the vdW/Casimir science and 
amount of published work warrant a separate review. 

We begin the discussion with the vdW regime (Sec. II). The most signihcant advances in the past decade have 
been in the development of novel Hrst principles methods for vdW calculations. Much of this progress has been 
motivated by the need for an accurate description of vdW interactions in materials as well as relevant experimental 
measurements. In the next section (Sec. Ill), we move on to larger separation scales and focus on emerging materials 
with Dirac-like spectra, such as graphene and systems with non-trivial topological phases. By summarizing results 
obtained via the Lifshitz theory, QED approach, and perturbative Coulomb interaction calculations we discuss how 
the Dirac spectra affect various characteristics of the vdW/Casimir force. The following two sections (Sec. IV and 
Sec. V) are devoted to Casimir interactions in structured materials. We discuss how the force can be manipulated 
via response properties engineering and non-trivial boundary conditions. For this purpose, we summarize not only 
important work in metamaterials, photonic crystals, and plasmonic nanostructures, but we also describe the progress 
in relevant computational tools. Biological materials are included in Sec. VI by highlighting results obtained via the 
Lifshitz and Hamaker theory calculations. Fluctuation phenomena for bio-systems are also discussed in light of other, 
Casimir-like phenomena. Much of this review is focused on the rapid expansion of theoretical and computational 
advances applied to the vdW/Casimir interactions. Although we concentrate on theoretical and computational work, 
key experiments giving us unprecedented insight into vdW and Casimir interactions are reviewed throughout the 
paper as well as in Sec VII. In the last section, we give our outlook for the future by discussing open problems in this 
field. 


II. AB INITIO METHODS FOR VAN DER WAALS FORCES 

Non-covalent interactions originating from correlated electron fluctuations between materials at separations on the A 
to a few nm scale play a key role in understanding their stability and organization. In recent years, important advances 
have been made towards computational methods for calculating vdW interactions with sufficient accuracy. These 
state-of-the-art methods are firmly based on a microscopic description of vdW interactions. Based on the treatment 
of the electron degrees of freedom of the atomistic system, we distinguish between two types of approaches: exact and 
approximate formulations of the many-body correlation energy. We discuss the essentials in terms of the adiabatic 
connection fluctuation-dissipation theorem (ACFDT), as both approaches rely on it. Based on the substantial evidence 
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accumulated over the last few years, we argue that ubiquitous many-body effects in the vdW energy are crucial for 
accurate modeling of realistic materials. The inclusion of these effects in first-principles calculations and comparative 
performance evaluation for a wide range of materials, including finite and periodic molecular systems, (hard) insulating 
and semiconducting solids, and interfaces between organic and inorganic systems are also discussed. 


A. Exact non-relativistic treatment of microscopic vdW interactions 


The exact energy of a microscopic system obtained via the solution of its Schrodinger equation seamlessly includes 
the vdW contribution. Explicitly solving the Schrodinger equation for more than a few electrons, however, is still 
a prohibitive task due to the complexity of the many-body problem. Therefore, first-principles modeling of realistic 
materials often starts with more tractable mean-field models, such as the Hartree-Fock approximation (HFA), or 
density-functional approximations (DFAs), which utilize the three-dimensional electron charge density, n(r), in lieu 
of the more complicated many-electron wavefunction. Unfortunately, these commonly utilized approximations are 
unable to describe the long-range electronic correlation energy and therefore fail to treat vdW interactions. 

The vdW energy is directly related to the electron correlation energy, Ec, which can be constructed exactly by 


invoking the ACFDT (Gunnarsson and Lundqvist, 1976 Langreth and Perdew, 1977) 


E. = -- 

2tt 


nOO nl 

/ duj dX Tr[{xx{r, r', iuj) - xo{r,r',iuj))v{r,r')], 
Jo Jo 


( 1 ) 


where ^^(r, r',iw) and are respectively the interacting and bare (non-interacting) response functions at 

Coulomb coupling strength A. Here, oj is the frequency of the electric field, u(r,r') = |r — r'|“^ is the Coulomb 
potential, and Tr denotes the spatial trace operator (six-dimensional integral) over the spatial electronic coordinates 
r and r'. The essential idea is that Eq. [^is an adiabatic connection between a reference non-interacting mean-field 
system with A = 0 and the fully interacting many-body system with A = 1 (Cunnarsson and Lundqvist 1976 Langreth 


and Perdew, 1977). The vdW contribution can be found from the so-obtained Ec in a tractable manner provided 


that a set of single-particle orbitals computed with DFAs or HFA can be used to construct xo(^,r^*w). This is still 
a formidable computational task for systems with thousands of electrons. In addition, approximations are needed to 
obtain XA(r,r',ia;) for 0 < A < 1. 

The power and significance of the ACFDT approach is that essentially all existing vdW modeling methods can 
be derived from approximations to Eq. For example, the widely employed pairwise approximation is obtained by 
truncating the ACFDT expression to second order in the perturbative expansion of the Coulomb interaction. The 
simple addition of inter-atomic vdW potentials that is used to compute the vdW energy in classical force fields and DFA 
calculations can be recovered from Eq. by further approximating the response function as a sum of independent 
dipole oscillators located at every nucleus in a given material (Tkatchenko et al. 2013). The vdW-DF approach 


originated by Langreth, Lundqvist, and collaborators (Cooper et al.. 2010 Dion et al. 2004 Lee et al. 2010) that has 


become widely used to correct semi-local DFAs can also be derived from Eq.jfjby making a local approximation to the 
response function in terms of the electron density and then employing second-order perturbation theory. However, 
the main shortcoming of all these rather efficient approximations is that they are unable to capture the non-trivial 
many-body effects contained in the interacting response function ^^(r, r',iw) as well as the infinite-order nature of 
the ACFDT expression in Eq. [l] 


B. Response functions and polarization waves 


The interacting response function is defined self-consistently via the Dyson-like equation 

XA = Xo + XoiXv + /r)XA, 


( 2 ) 


which contains the exchange correlation kernel /“'^(r, r', iw), an unknown quantity which must be approximated in 
practice. Neglecting the explicit dependence of on the coupling constant allows for an analytic integration over A 
in Eq. |p^, and forms the basis for the most widely employed approximation, namely the random-phase approximation 


(RPA) (Bohm and Pines 1953 Cell-Mann and Brueckner, 1957). 


The non-interacting response function can be obtained using the Adler-Wiser formalism ( Adler[ 1962[ Wiser, 1963), 
given a set of occupied and unoccupied electronic orbitals with corresponding energies {e^} and occupation 
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numbers {fi} determined from semi-local DFT, Hartree-Fock, or hybrid self-consistent field calculations, be., 

\Ji Jj. 

^3 


Xo{v,r',iuj) = J2{fi - fj)- 


ei - e,- 


• ILU 


( 3 ) 


This mean-field xo can exhibit relatively long-range fluctuations (polarization waves), the extent of which is determined 
by the overlap between occupied and rather delocalized unoccupied electronic states used in Eq. In this framework, 
the fluctuations in xi may be shorter-ranged than in xo^ especially in 3D solids where the Coulomb interaction leads 
to significant screening effects. The situation is generally very different in anisotropic nanostructured materials, where 
the Coulomb interaction might lead to so-called anti-screening effects, i.e. significantly farsighted polarization waves 
(see Fig. [^a) for illustration). 

So far, the general understanding of polarization waves comes from coarse-grained approximations to the density- 
density response function. For example, Dobson et al. (Dobson et al. 20061 found that the asymptotic vdW 


interaction between two low-dimensional metallic objects differs qualitatively from the commonly employed sum-over- 
pairs expressions. Another example is the vdW graphene-graphene interaction energy decays as d~^, instead of the 
conventionally expected dr'^ power law. Here, however, many-body renormalization of the Dirac graphene carriers 
beyond the RPA might lead to vdW interaction power law between d~^ and d~^ (Dobson et al. 2014) (also discussed 
in Sec. IIIC). This is a matter of ongoing debate. 

For some time it was assumed that complete delocalization of fluctuations is required to identify interesting de¬ 
viations from the otherwise pairwise-additive behavior. However, Misquitta et al. (Misquitta et al.. 2014 2010) 


demonstrated that semiconducting wires also exhibit unusual asymptotics, which becomes more pronounced with the 
decrease of the band gap ( Misquitta et al.\ 2014 2010| . In this case, the vdW interaction exhibited a power law of d~^ 
at large but finite distances, converging to the standard d~^ behavior for large inter-wire separations. Ambrosetti et 
al. (Ambrosetti et al. 2014b) have analyzed the spatial extent of dipole polarization waves in a wide range of systems 


and demonstrated a continuous variation of the power law for finite distances between ID wires and 2D layers with 
visibly enhanced non-local responses due to the collective many-body effects. Such relative farsightedness of vdW 
interactions provides an avenue for appropriately tuning the interactions between complex polarizable nanostructures. 

Another way to understand polarization waves in materials consists in studying the renormalization (non-additivity) 
of polarizability and vdW coefficients for different systems as a function of their size and topology. Ruzsinzsky et al. 


(Ruzsinsky et al.. 2012) modeled the polarizability of fullerenes employing a hollow shell model with a hnite thickness. 


They demonstrated that the polarizability scales superlinearly as a function of fullerene size. This leads to a super¬ 
quadratic increase in the vdW Cg coefficients, clearly demonstrating the importance of long-range fluctuations. Tao 
and co-workers extended these findings to a wide range of nanoclusters (Tao and Perdew, 2014). Recently, Gobre 
and Tkatchenko (Gobre and Tkatchenko 2013) studied the dependence of carbon-carbon vdW coefficients for a 


variety of carbon nanomaterials and they found that vdW Cq coefficients could change from 20 hartree-bohr® to 
150 hartree-bohr® depending on the dimensionality, topology, and size of the carbon nanostructure. This clearly 
demonstrates the extreme non-additivity of vdW interactions in low-dimensional materials, and highlights the need 
to include collective effects in vdW interactions when modeling the self-assembly of such nanostructures. 


C. Approximate microscopic methods for van der Waals interactions 


Since the vdW energy is a tiny part of the total energy of a many-electron system, vdW methods have to be 
coupled to an underlying electronic structure method that provides an adequate treatment of hybridization, charge 
transfer, electrostatics, and induced polarization, among other electronic structure effects. Density-functional theory 
(DFT) with approximate exchange-correlation functionals provides an optimal approach in this regard. DFT is able to 
correctly describe short-range quantum-mechanical interactions, and also treats classical electrostatic and polarization 
effects rather accurately. The total energy Et of a many-electron system is 


Et — Akin + Aes + Ax -|- Ac 


( 4 ) 


where Akin is the electronic kinetic energy (corresponding to mean-field kinetic energy in the Kohn-Sham framework), 
Aes is the electrostatic energy (including nuclear repulsion, electron-nucleus attraction, and Hartree electron-electron 
repulsion), and Ax and Ac are the non-classical exchange and correlation terms, respectively. Most DFT methods 
utilize semi-local approximations by using information about the electron density (local density approximation, LDA) 
and its gradients (generalized gradient approximation, GGA). Other approaches are based on the Laplacian of the 
electron density in the so-called meta-GGA functionals (Sun et al. 2013 Zhao and Truhlar 2008). It may also 
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(c) Two-point functionals 


(d) Fragment-based approaches 



FIG. 2 (Color online) Schematic representation of first principles methods for the (a) exact formulation of the electronic 
correlation energy Ec from the adiabatic connection fluctuation-dissipation theorem; (b) formulation based on coupled dipolar 
fluctuations, such as the many-body dispersion (MBD) methods ( [Ambrosetti et aT] |2014^ [Tkatchenko et aZ.[ |2012b[ ); (c) Ec 
using two-point functionals obtained by approximating the non-homogeneous system with a homogeneous-like response; and 
(d) fragment-based correlation energy obtained from multipolar expansions. 


be advantageous to include a certain amount of exact Hartree-Fock exchange in DFA, leading to so-called hybrid 
functionals. 

We note that from ACFDT and Ec are non-local (Eq. [^. The correlation energy Ec, which is of relevance to 
the vdW interaction, can be written as Ec = E^\ + E^y, where Es\ is the semi-local correlation energy and E^y is the 
non-local part. The fact that such partition is not unique has led to a flurry of heuristic approaches that aim to 
construct a reliable approximation to the full electronic correlation energy. The different classes of methods for the 
non-local correlation energy are schematically shown in Fig. and are summarized in what follows. 


1. Two-point density functionals for vdW interactions 


Obtaining an exact expression for XA(r, r', w) in general is not possible. However, for a 3D homogeneous electron gas 
the correlation energy can be written exactly in terms of the electron density n(r). Approximating the polarization of 


a non-homogeneous system assuming homogeneous-like response is possible in certain situations (Dobson and Dinte 
1996[ [Rapcewicz and Ashcroft 1991). These ideas have led to the derivation of the vdW-DF approach. In addition 


to the approximation of the interacting polarizability as a local quantity, one also takes a second-order approximation 
in Eq. [^assuming x\ = Xi- Thus the non-local correlation energy is obtained as 


e: 


vdW-DF 

nl 


= J n{r)K{r,r')n{r')drdr' 


( 5 ) 


where K{r,r') is a “vdW propagator” (Fig. [^c)). Note that Eq. constitutes a great simplification over the exact 
Eq.[T] since only n(r) and its gradient (utilized in K) are required. 

The original implementation of this additive non-local correlat ion energy to the total DFT energy was proposed to 


couple t;™" to a revised Perdew-Burke-Ernzerhof functional (Zhang and Yang 19981, the rationale being that this 


functional yields repulsive binding-energy for prototypical vdW-bound systems, such as rare-gas dimers (Dion et al. 


2004). A follow-up implementation using the DFT functional PW86 (Lee et al.. 2010) has generated a revised, vdW- 


DF2 functional. While the vdW-DF2 approach was shown to perform much better for intermolecular interactions. 


its behavior at vdW distances is significantly deteriorated when compared to vdW-DF (Vydrov and T. Van Voorhis 
2010). Specifically, while the vdW Cg coefficients in the vdW-DF method are accurate to 19%, the error increases to 


60% when using vdW-DF2. These approaches illustrate the challenging problem of balancing between semi-local and 
non-local interactions in a meaningful manner. 

Following the success of the vdW-DF approach, Vydrov and Van Voorhis (VV) provided a significantly simplified 
vdW functional derivation and revised the definition of local polarizability, by employing a semiconductor-like dielectric 
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function along with the Clausius-Mossotti relation between polarizability and dielectric function (Vydrov and T. Van 


Voorhis 2009 2012). The VV approach requires one parameter for the local polarizability and a second one for the 


coupling between the non-local vdW energy with the parent DFA approach. The VV functional was assessed with a 
wide range of semi-local and hybrid functionals, and by tuning these two parameters, it yielded remarkable performance 


for intermolecular interactions compared to benchmark data from high-level quantum-chemical calculations (Vydrov 


and T. Van Voorhis, 2012). Other approaches, such as the C09x functional of Cooper (Cooper 2010) and the “opt” 
family by Klimes and Michaelides (Klimes et al. 2010 2011), rely on the same definition in Eq. however the coupling 
with the DFA is revised by adjusting one or more parameters in the semi-local functional. The ’’opt” functionals 
parameters in particular were adjusted to a benchmark database of intermolecular interaction energies showing a good 
performance for cohesive properties of solids (Klimes et al. 2011). There are indications, however, that the ’’opt” 


functionals overestimate the binding in larger and more complex molecular systems (Klimes and Michaelides 2012). 


These recent developments have led to many novel insights into the nature of vdW interactions. However, the 
drastic approximations in tennis of the additive polarizability and the dependence on the electron density 

on two points only must be assessed carefully for realistic materials. The neglected non-additive effects can play a very 


important role in many systems (Ambrosetti et al., 2014b Dobson et al., 2006 Gobre and Tkatchenko 2013 Misquitta 


et al. 2014 2010 Ruzsinsky et ol. 2012 Tao and Perdew, 2014 Tkatchenko 2015). Also, the neglected three-body 


2012 

Donchev 

et al. 

2012b 

). 7 


2014a 


DiStasio, Jr. et al. 


Donchev 2012[ |Kronik and Tkatchenko 2014 Marom et al. 20131 [Shtogun and Woods[ 2010[ Tkatchenko 


without substantially increasing their cost. One possibility is going towards RPA-like approaches, but this would 
mean departing from a pure density functional picture. Another possibility entails further coarse graining of the 
system to a fragment-based description, the progress of which is summarized below. 


2. Fragment-based metliods for vdW interactions 


Fragment-based methods can be traced back to the original work of London ( London[ |1930 ), in which case utilizing 
second-order perturbation theory for the Coulomb interaction the dispersion energy between two spherical atoms A and 
B can be obtained. The London expression, often using just the dipolar term Cq^ab/R\b^ i® the basis for calculating 


vdW dispersion energies in a wide range of atomistic methods, including Hartree-Fock calculations (Ahlrichs et al. 
1977 Hepburn and Scoles 1975), DFA calculations (Grimme 2006 Grimme et al. l IMol [Johnson and Becke[ |2005 


et al. 


Steinmann and Gorminboeuf 2011 Tkatchenko and Scheffler] 2009) and quantum chemistry methods (Tkatchenko 


l‘2\> ' 

2009). is valid only at large separations for which the overlap between orbitals of atoms A and B can be 


neglected. At shorter sepations the overlap naturally reduces the interaction ( 

Koide 

1976 

, which can be conviniently 

included by a damping function 

' Grimme 

2006 

Grimme et al. 

2010 

Johnson and Becke 

2005 

Tang and Toennies 

1984 

Tkatchenko and Scheffler 

2009) 


p(2) _ 

-^vdW “ 


6 , 8 , 10 ,. 


fd.n{BAB ^ Bc,Ab'} 


Cr, 


,AB 


R 


( 6 ) 


AB 


where fd,n{RAB,Rc.AB) is the damping function that depends on a cutoff radius Rc,ab (Fig- [^d)). This type of 
approach can be quite useful in DFA-GGA functionals, such as PBE (Perdew et al. 1996), which perform very well 
for chemical bonds. In this case, the dipolar approximation to Eq. |6|is sufficient (Elstner and Hob'za||2001 Grimme 


to 

o 

o 

Wu and Yang 

2002 

Wu et al. 

2001 

Zimmerli et al. 


tremendous developments. In particular, Grimme (Grimme 2006) published a set of empirical parameters for a range 


|‘ 


2004|. These ”DFT-D” approaches have experienced 


of elements and demonstrated that the addition of dispersion energy to a wide range of functionals yields remarkably 
accurate results for intermolecular interactions. In the latest DFT-D3 method, Grimme has extended his empirical 


set of parameters to cover elements from H to Pu (Grimme et al ., 2010). 


Gonsiderable efforts have been dedicated towards determining vdW parameters directly from electronic structure 
calculations also. Becke and Johnson (Johnson and Becke, 2005) developed an approach based on the exchange- 


hole dipole moment (XDM) to determine vdW coefficients, which can be computed by using Hartree-Fock orbitals. 
Later, Steinmann and Gorminboeuf presented an alternative derivation based on electron density and its first and 
second derivatives ([Steinmann and Gorminboeuf 2011). The alternative derivation of the fragment-based vdW- 


DF functional by Sato and Nakai has demonstrated an interesting connection (and potential equivalence) between 
fragment-based methods and explicit non-local functionals (Sato and Nakai 2009 2010). Tkatchenko and Scheffler 


(TS) developed a DFA based approach to determine both Cq^ab coefficients and Rc,ab radii as functionals of the 













































































































































































































































electron density (Tkatchenko and Scheffler 20091, which implies that the vdW parameters respond to changes due 


to hybridization, static charge transfer, and other electron redistribution processes. The TS approach demonstrated 
that by utilizing the electron density of a molecule and high-level reference data for the free atoms, it is possible to 
obtain asymptotic vdW coefficients with accuracy of 5.5%, improving by a factor of 4-5 on other existing approaches 
at the time. Bucko and co-workers have pointed out that an iterative Hirshfeld partitioning scheme for the electron 
density can significantly extend the applicability of the TS method to ionic materials (Bucko et al. 2014, 2013a|. 

This field is still developing at a rather quick pace, therefore revised and completely new fragment-based approaches 
are still being introduced. 


3. Efficiently beyond pairwise additivity; Explicit many-body vdW methods 


In more complex and heterogeneous systems, it is necessary to go beyond the simple additive models and further 
efforts of atomistic vdW modeling must be directed towards the inclusion of many-body effects. In principle, RPA 
using DFA orbitals provides a good model, however the dependence of xo on the exchange-correlation functional and 
the high computational cost in the xo computations may be limiting factors. The main challenge is to construct reliable 
approximations for the long-ranged vdW correlations, since the short-ranged correlation effects are well accounted for 
in DFA. Therefore, the full xo(^,r^^a;) is often unnecessary as is the case in nonmetallic or weakly metallic systems. 
In such situations, it is possible to describe xo by a set of localized atomic response functions (ARFs), which can be 
constructed to accurately capture the electronic response beyond a certain cutoff distance. 


2009 

Donchev 

2012 Liu et al. 

2011b Shtogun and Woods 

2010 

I, only recently has this idea been extended to non- 

local vdW interactions in realistic materials ( 

DiStasio, Jr et al.\ 

2014 

). For this purpose, spatially-extended ARFs 


that increase the applicability of the model to include close contact have been used within the so-called many-body 
dispersion (MBD) method, schematically illustrated in Fig. [^b) ( [Ambrosetti et al. 2014c Tkatchenko et al. 2012bI. 
Within this approach, each p-th atom in the material is represented by a single dipole oscillator with a frequency- 

O' 

dependent polarizability ap{iuj) = ; where app is the static polarizability and lUp^ is an effective excitation 


(or resonant) frequency. The bare ARF response then is written as 




, r', iuj) = —ap{iu!)V {r — Rp) 0 Vr'5^(r' — Rp), 


( 7 ) 


where Rp is the location of the p-th atom and 0 signifies a tensor product. The bare response function for a collection 
of atoms follows simply as the direct sum over the individual ARFs, Xo(r, r', iuj) = Xo,p(r, 0 Xo.( 3 (r, r', iuj) 0 • • •. 

The ARF response contains the infinite-order correlations from the start and it can be used in Eq. to calculate 
the interaction energy. It has been demonstrated that the RPA correlation energy is equivalent to the exact diag- 


onalization of the Hamiltonian corresponding to ARFs coupled by a long-range dipole potential (Tkatchenko et al. 


20131. Using second-order perturbation theory one also recovers the well-known pairwise-additive formula for the 


vdW energy (Tkatchenko et al. 2013). 

We note that the solution of Eq. for a model system of AREs yields an expression for the long-range correlation 
energy beyond what would simply be called “vdW dispersion energy” in the traditional London picture (London, 19301. 


Even for two atoms described by dipole-coupled ARFs, the correlation energy contains an infinite numbers of terms 
Cn,AB/R^ ab- The polarizability of the combined AB system in general is not equal to the sum of polarizabilities 
of isolated atoms A and R, and higher-order correlation terms account precisely for this fact. It has been found 
that the convergence of the perturbative series expansion in Eq. can be extremely slow, especially for systems 
which have either high polarizability density or low dimensionality. This is clearly illustrated by the binding energy in 
supramolecular complexes or double-walled nanotubes, for which even 8-body terms make a non-negligible contribution 


to the correlation energy on the order of 2-3% (Ambrosetti et al.. 2014a I. 


D. Applications of atomistic vdW methods to materials 


Our discussions above show that the developments of novel many-body methods and understanding of many-body 


effects in the vdW energy is an area of significant current interest ( 

Ambrosetti et al. 2014c 

Modrzejewski et al. 

2014 

Otero-de-la Roza and Johnson 

2013 

Silvestrelli 2013 

Tkatchenko et a/.||2012b 

|. This is 

highly motivated from an 


experimental point of view as well. Being able to describe vdW interactions in different systems is highly desirable 
in order to explain existing and predict new experimental findings. Figsummarizes typical results from ab initio 
calculations, as described below, for cohesion energies and error ranges as compared to available reference data. 
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Gas-phase 

molecules 

10-20 mcV(5-I0%) 

Pristine 

solids 



0.2-0.3 eV/atom 


Supramolecular 

systems 



Molecular 

crystals 

0.04-0.1 eV(5-15%) 

Organic/inorganic 

systems 


Defects 
in solids 



up to 0.4 eV - I eV per benzene ring 


FIG. 3 (Color online) Various types of materials with calculated cohesion energies and errors as compared to available reference 
data. The top row of valnes shows typical errors of atomistic vdW methods compared to benchmark bending energies, while 
the bottom row of values shows contributions of the vdW energy to the binding energy of the corresponding materials. 


1. Finite and periodic molecular systems 


For smaller molecules, high-level quantum-chemical benchmarks using coupled-cluster calculations are now widespread 
urecka et al. 2006 Rezac et al. 2011[ Takatani et al. 20101. In particular, the coupled-cluster method with single, 
double, and perturbative triple excitations [CCSD(T)] is currently considered as the “gold standard” of quantum 
chemistry. For yet larger molecules (up to 200 light atoms), it is possible to carry out Diffusion Quantum Monte 
Carlo (DQMC) calculations ( Ambrosetti et alj [2014^ Benali et al.. 2014[ ) using massively-parallel computer ar¬ 
chitectures. DQMC calculations in principle yield the exact solution (within statistical sampling accuracy) for the 
Schrodinger equation within the fixed-node approximation (Foulkes et all 12001 1. Examples of small molecules bench¬ 


mark databases are those for the S22 (Jurecka et al.. 2006 Takatani et ail 20101 and S66 (Rezac et all 20111 dimers, 


containing 22 and 66 dimers, respectively. For supramolecular systems, the S12L database has been recently intro- 


duced by Grimme ( 

Grimme 

20121 and benchmark binding energies for 6 of these 12 complexes have been calculated 

using DQMC ( 

Ambrosetti et al. 

2014a 

). For extended periodic molecular crystals, one can rely on experimental 


lattice enthalpies, extrapolated to 0 K and with zero-point energy subtracted. Two databases, C21 (Otero-de-la Roza 


and Johnson 2012) and X23 (Reilly and Tkatchenko 2013aI, have been recently introduced for molecular crystals. 


Both of these databases include molecular crystals bound primarily by either hydrogen bonds or vdW dispersion, 
including a few crystals with mixed bonding nature. The X23 database extended the C21 one and improved the 
calculation of vibrational contributions required to convert between experimental sublimation enthalpies and lattice 
energies. 

Initially, the development of atomistic methods for vdW interactions has been largely driven by their performance 
for small molecules in the S22 and S66 databases. Currently, the best methods are able to achieve accuracies of 10-20 
meV (better than 10%) in the binding energies compared to reference CCSD(T) values (Fig. [^. The errors are 
due to inaccuracy in the asymptotic vdW coefficients, empirical parameters in damping functions, and errors in the 
exchange-correlation functional. 

Due to such rather uniform performance of different methods for small molecules, the focus has shifted to assessing 
the performance for larger systems. Here, in fact, the differences are more prominent, because the vdW energy makes a 
much larger relative contribution to cohesion. For example, for polarizable supramolecular systems, such as the “buck- 


yball catcher” complex, pairwise dispersion corrections overestimate the binding energy by 0.4-0.6 eV (Tkatchenko 


et al.\ 2012a) compared to reference DQMC values. Only upon accurately including many-body dispersion effects 
one obtains results within 0.1 eV from the best available benchmark (Ambrosetti et al. 2014aI. So far, vdW-DF 


functionals have not been applied to study binding energetics in the S12L database. 

For periodic molecular crystals, some pairwise and many-body fragment-based methods are able to achieve remark¬ 


able accuracy of 40-50 meV per molecule (5% mean absolute relative error), compared to experimental results (Reilly 
and Tkatchenko 2013a|b ). Since the difference in lattice energies between various available experiments is on the same 


order of magnitude, this highlights the mature state of vdW dispersion corrections to DFA. The vdW-DF2 approach 
yields a somewhat larger error of « 70 meV (7.5%) on the C21 database (Otero-de-la Roza and Johnson 20121. 
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Understanding the performance of different vdW-inclusive methods for large molecular systems is still a subject of 
ongoing research. Some of the pairwise correction approaches have been specifically fitted to periodic systems, trying 
to mimic many-body screening effects by changing the short-range damping function. This procedure seems to work 
well for certain molecular crystals with high symmetry, but this is obviously not a transferable approach. 

Many-body vdW correlations become even more relevant for the relative energetics of molecular systems, which are 
essential to predict the correct polymorphic behavior of molecular crystals. Marom et al. have demonstrated that only 
upon including many-body effects one is able to correctly reproduce the structures and relative stabilities of glycine, 
oxalic acid and tetrolic acid (Marom et al. 2013). Another interesting example is the aspirin crystal, for which a 


long-standing controversy has been about the relative stability of polymorphs form I and form II ([Ouvrard and Price 


2004). Reilly and Tkatchenko have recently demonstrated that the stability of the most abudant form I arises from 


an unexpected coupling between collective vibrational and electronic degrees of freedom (dynamic plasmon-phonon 
coupling) (Reilly and Tkatchenko 2014). In this case, many-body vdW correlations renormalize phonon frequencies 
leading to low-frequency phonon modes that increase the entropy and ultimately determine the stability of this 
ubiquitous form of aspirin in comparison to the metastable form II. Furthermore, the bulk and shear moduli of both 
forms are substantially modified and become in better agreement with experiments when calculated with DFA-I-MBD. 
The aspirin example illustrates how the inclusion of many-body vdW effects may lead to novel qualitative predictions 
for the polymorphism and elastic response of molecular materials. 


2. Condensed materials 


For hard solids (ionic solids, semiconductors, and metals), the role of vdW interactions was considered to be 


negligible for a long time, as judged for example by classical condensed-matter textbooks (Ashcroft and Mermin 


1976 Kittel 1986). The rather strong cohesion in hard solids stems from covalent and metallic bonds, or from 


classical Coulomb interaction between localized charges (Fig. [^. Early estimates of vdW interactions in hard solids 


varied substantially from being negligible to being very important (Ashcroft and Mermin 

1976 

Mayer 

1933 

Rehr 

et al. 

1975 

Richardson and Mahantyf 1977 Tao et al.. 2010|). Recently, this issue has been systematically revisited 


by employing DFA with vdW interactions. Zhang et al. ( [Zhang et al.\ |2011| ) demonstrated that long-range vdW 
interactions account for « 0.2 eV/atom in the cohesive energy for Si, Ge, GaAs, NaCl, and MgO, and 9-16 GPa in 
the bulk modulus. This amounts to a contribution of 10-15% in the cohesive energy and bulk modulus - far from 
being negligible if one aims at an accurate description of these properties. Klimes and co-workers applied their “opt” 


functionals based on the vdW-DF approach to a large database of solids (Klimes et al. 2011) finding that vdW 


interactions play an important role for an accurate description of cohesive properties. Overall, their conclusion is that 
vdW interactions allow improving the performance of many different xc functionals, achieving good performance for 
both molecules and hard solids. 

Because vdW interactions are important for absolute cohesive properties of solids, any property that depends on 
energy differenees is also likely to be influenced by vdW effects. Therefore, vdW interactions will often play an 
important role in the relative stabilities of different solid phases, phase transition pressures, and phase diagrams as 


Moellmann et al. 

2012 

1, ice ( 

Santra et al. 2011 


phases of the oxidized Cu(llO) surface (Bamidele et al. 2013), and alkali borohydrites (Huan et al. 2013). 


2011), different reconstructed 


The properties of many solids are substantially affected by the presence of simple and complex defects, such as 


neutral and charged interstitials and vacancies (Freysoldt et al. 2014). The formation of defects entails a modification 


of polarization around defect sites and this can have a substantial effect on the contribution of vdW energy to the 
stability and mobility of defects. Gao and Tkatchenko have demonstrated that the inclusion of many-body vdW 
interactions in DFA improves the description of defect formation energies, significantly changes the barrier geometries 
for defect diffusion, and brings migration barrier heights into close agreement with experimental values (Gao and 


Tkatchenko 2013). In the case of Si, the vdW energy substantially decreases the migration barriers of interstitials 


and impurities by up to 0.4 eV, qualitatively changing the diffusion mechanism (Gao and Tkatchenko 2013). Recently, 


the proposed mechanism has been confirmed by explicit RPA calculations (Kaltak et al. 2014). 


3. Interfaces between molecules and solids 


The predictive modeling and understanding of hybrid systems formed between molecules and solids are an essential 
prerequisite for tuning their electronic properties and functions. The vdW interactions often make a substantial 
contribution to the stability of molecules on solids (Tkatchenko et al.. 2010). Indeed, until recent developments for 
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efficiently incorporating the long-range vdW energy within state-of-the-art DFAs, it was not possible to study the 
structure and stability of realistic interfaces (Liu et al.. 2014). Exposed surfaces of solid materials are characterized by 
formation of collective electronic states, thus the long-range screening effects should be treated at least in an effective 
way, as done for example in the DFT-|-vdW®“''^ method that accounts for the collective electronic response effects by a 
combination between an interatomic dispersion expression and the Lifshitz-Zaremba-Kohn theory (Ruiz et al. 2012|). 


This method was demonstrated to be reliable for the structure and stability of a broad class of organic molecules 
adsorbed on metal surfaces, including benzene, naphthalene, anthracene, diindenoperylene, Ceo, and sulfur/oxygen- 
containing molecules (thiophene, NTCDA, and PTCDA) on close-packed and stepped metal surfaces, leading to an 
overall accuracy of 0.1 A in adsorption heights and 0 . 1 - 0.2 eV in binding energies with respect to state-of-the-art 
experiments (Liu et al. 2014|). 


A particularly remarkable finding is that vdW interactions can contribute more to the binding of strongly bound 


molecules on transition-metal surfaces than they do for molecules physisorbed on coinage metals (Carrasco et al. 


2012 Liu et al. 2012). The accurate inclusion of vdW interactions also signihcantly improves molecular tilting angles 


and adsorption heights, and can qualitatively change the potential-energy surface for adsorbed molecules with flexible 


functional groups. Activation barriers for molecular switches (Liu et al. 2013) and reaction precursors (Rodriguez- 


Reyes et al ., 2014) are modified as well. Ongoing work concentrates on understanding the interplay between many-body 


effects within the solid material and collective effects within the adsorbed molecular layers. 


DIRAC MATERIALS BEYOND ATOMIC SCALE SEPARATIONS 


Without a doubt, novel ab initio methods have advanced our understanding of vdW interactions between systems 
at atomic scale separations. When the objects are taken further apart, however, other approaches become more ap¬ 
propriate. Dispersion interactions involving objects with macroscopic dimensions at distances for which the electronic 


distribution overlapping effects are not important are typically described by the Lifshitz formalism (Dzyaloshinskii 


et al. 1961 Lifshitz 1956), which has been the mainstream theory for conventional metals and dielectrics with pla¬ 


nar extensions for over several decades. New materials with Dirac spectra are emerging, however, and the Lifshitz 
approach is an excellent tool to capture the signatures of the Dirac carriers in the vdW/Casimir interactions. 


A. Lifshitz formalism 


A generalized Lifshitz formula can be obtained from the ACFDT expression in Eq. for distance separations 
larger than several nm-s, where the overlap of the electronic distribution residing on each object can be neglected. In 
this case, the response properties are independent of each other, thus they are described by the individual response 
functions Xq ’ ' (1,2 denote the two objects) and the mutual Coulomb potential can be taken as a perturbation. 
Although Xo ’ do not include electronic correlations from the overlap, they contain the electronic correlations within 
each object. When Xo ’ are calculated via the RPA approach and the mutual Coulomb interaction is described by 
the RPA ring diagrams (Dzyaloshinskii e< aZ. 1961 Fetter and Walecka 1971 Lifshitz 1956), the ACFDT expression 
(Eqp]) is transformed to the non-retarded Lifshitz formula given in Fourier basis for planar homogeneous objects 


(Despoja et al. 2007 Dobson and Gould 2012) 


= -n 


dkji 

Wf 


doj 


ln[l - Xo^ (k||, iu})Vi2 (ky , iuj)xo^’ ^ (ky , iuj)] 


( 2 )/ 


( 8 ) 


where ky is the 2D wave vector and Vi 2 (ky,ia;) is the Coulomb interaction between the two objects. 

Distance separations on the order of sub-/rm and /r m scales are characteristic for the Casimir regime, where 
retardation becomes prominent. For such separations one has to include all photon interactions being exchanged with 
the finite speed of light c. The Casimir energy can be derived using scattering methods by solving the boundary 
conditions arising from the electromagnetic Maxwell equations. The interaction energy can written in the form 


Lambrecht et al. 

2006 

Lambrecht and Marachevsky 

2008 Rahi et al. 

2009 


=—h [^ [ I InDetp — Mi(ia;)K 2 (ia;)e 




-+kS 


(9) 


where Ki _2 are the reflection matrices of the individual objects evaluated at imaginary frequencies. The reflection 
matrices describe appropriate boundary conditions and they are expressed in terms of the macroscopic response 
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properties of the objects. Eq. can be obtained equivalently via QED techniques relying on the evaluation of the 
Maxwell stress tensor whose components represent the vacuum expectation of the electromagnetic field and they 


are given in terms of the dyadic Green’s function (more details on this approach are found in Sec. V) (Buhmann 


and Welsch 

2007 

Dzyaloshinskii et al. 

1961 


Utilizing the fluctuation-dissipation theorem and standard complex 
contour integration techniques, Matsubara frequencies iuin = in2'KkBT/H are introduced in the description. As 
a result, the Casimir interaction energy can be cast into a temperature-dependent form using the relation 

h —>• prime in the sum means that the n = 0 term is multiplied by 1/2). We further note that 

setting c = oo in Eq. ^the non-retarded Lifshitz expression (Eq. is recovered. 

Being able to utilize independently calculated response properties with different models (including RPA or the 
Kubo formalism) or even use experimental data in Eqs. 8[9 has been especially useful for the versatility of the Lifshitz 
macroscopic approach. In addition, the Matsubara frequencies give the means to take into account temperature in the 
interaction unlike the ab initio methods (Sec. II), which calculate the vdW energy at zero temperature. Much of the 
progress in theoretical and experimental work concerning typical metals and dielectrics interactions, captured by the 


macroscopic vdW/Casimir approach, has been summarized in several books and recent reviews (|Bordag et al. 2009b 
Buhmann and Welsch[ 2007 Dalvit et al. 2011a Klimchitskaya et al. 2009 Lamoreaux 2005[ j^rsegian 20061. 

The materials library is expanding, however. A new subset of systems, characterized by Dirac fermions in their low- 


energy spectra, has emerged recently (Welding et al. 20141. This distinct class of materials has properties markedly 


different from the ones of conventional metals and semiconductors whose fermions obey the Schrodinger’s equation. 
Recent discoveries have shown that there are many types of systems with Dirac nodes in the band structure, including 
2D graphene, topological insulators (TIs) and Weyl semimetals. Research efforts on vdW/Casimir interactions involv¬ 
ing graphene and related systems have shown that fluctuation-induced phenomena are strongly influenced by the Dirac 
nature of the carriers. As discussed previously, ab initio calculations have been indispensable for the demonstration of 
atomic registry-dependent effects, unusual scaling laws, farsightedness and many-body nature of their vdW interaction 

2013a|b 


in graphitic systems (Bucko et al. 2014 


Gobre and Tkatchenko 2013 Lebegue et al. 2010 Shtogun and 


Woods 2010). Nevertheless, it is very important from a fundamental point of view to consider regimes where the 
dispersion interactions are determined primarily by the low-energy Dirac carriers. Unlike ab initio methods which 
take into account the entire band-structure of the interaction materials, the Lifshitz/Casimir formalism relying on 
response functions calculated via low-energy models gives us an excellent opportunity to study the emergent physics 
of the Dirac carriers in vdW/Casimir forces. 


B. Basic properties of graphene nanostructures 


2004), significant progress has been made towards understand- 


After the discovery of graphene (Novoselov et al. 
ing its properties. For example, basic science in terms of the 2D Dirac-like nature (Neto et al. 


transport (Das Sarma et al.\ 2011), collective effects due to electron-electron interactions (Kotov et al. 


2009), electronic 
|2012[ ), and 


spectroscopy (Basov et al.\ 2014) have been studied. Quasi-ID allotropes, such as carbon nanotubes (GNTs) and 


graphene nanoribbons (GNRs) are also available ( 

Ma et al. 

2013 

Saito et al. 

1998 

), and key scientific breakthroughs 

have been summarized (Gharlier et al. 

2007 

Yaziev 201( 

). In addition to the internal properties, understanding 


primary importance. Much progress in the past several years has been achieved towards learning how such dispersive 
forces are influenced by the graphitic internal properties and external factors, such as temperature, doping, and ap¬ 
plied fields. This knowledge is relevant for a variety of phenomena including the formation and stability of materials 
and composites, adsorption, manipulation of atoms, and operation of devices, among others. Since the description of 
the interaction via the macroscopic Lifshitz/Gasimir formalism depends upon the low-energy electronic structure and 
optical properties, here we provide an overview of the relevant characteristics of graphene, GNTs, and GNRs. 

Graphene is a 2D atomic layer composed of hexagonally oriented rings (Fig. |^. Many of its properties can be 
captured by a nearest-neighbor tight-binding model within the first Brillouin zone with two inequivalent K points 
at K = (±^^^,0) (a is the graphene lattice constant), which describes a tt valence bonding band with one electron 
and an empty tt* anti-bonding band with one hole. The linearization of the energy spectrum around the K centered 
valleys yields the low-energy massless chiral Dirac-like Hamiltonian in 2D 


Hgr = hvFcr • q — (10) 

where /r is the chemical potential and a is the 2D spinor. The nonzero spinor components ay are the Pauli matrices, 
which refer to the graphene pseudospin rather than the real spin. The energy spectrum Egr is linear with respect to 
the wave vector q = k — K according to Egr = hvpq with the electronic group velocity being vp = V^tQa/2h ~ 10® 
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FIG. 4 (Color Online) (a) A graphene layer is an atomically thin sheet of honeycomb carbon atoms. Zigzag (green) and armchair 
(blue) graphene nanoribbons can be realized by cutting along the specified edges. Carbon nanotubes can be obtained by folding 
along the chirality vector Ch, determined by the indices n and m and the lattice unit vectors 01 , 2 - The chirality index (n,m) 
uniquely specifies each nanotube with two achiral examples shown: (b) zigzag (m,0) and (c) armchair (n,n). Alternatively, 
folding a zigzag nanoribbon along the axial direction results in an armchaired nanotube, while folding an armchair nanoribbon 
gives a zigzag nanotube. 


m/s {to is the nearest-neighbor tight-binding hopping integral). It is interesting to note that although graphene was 
synthesized not long ago, theoretical insight in terms of the low-energy massless Dirac-like Hgr was discussed much 


earlier by (Semenoff 1984), who expanded upon the tight-binding description introduced by (Wallace 1947). 


The tight-binding model can be extended to CNTs as well. Imposing periodic boundary conditions around the 
cylindrical circumference, the corresponding wave functions are zone-folded, meaning that the wave vector in the 
azimuthal direction takes a set of discrete values. The nomenclature of CNTs is described via a chirality vector 
Ch = nsLi + ma 2 , where n,TO are integers. As a result, single-walled CNTs are denoted via a chirality index {n,m) 
with the achiral nanotubes labeled as armchaired (n, n) or zigzag (n, 0), as shown in Fig. The CNT energy 


bands can also be found (Mintmire et al. 1992 Tasaki et al. 1998) with zone-folding boundary conditions leading 


to a chirality dependent energy spectrum. The tight-binding energy band structure for GNRs, on the other hand, is 
obtained by requiring the wave function be periodic along the CNR axis and vanish at the edges, which introduces 


edge dependent (zigzag or armchair) phase factors in the energy spectra (Akhmerov and Beenakker 2008 Brey and 


Fertig 2006 Sasaki et a/., 20111. 


The electronic structure of graphene systems determines their optical response properties - key components for 
the vdW/Casimir calculations. The Dirac-like carriers have profound effects on how electromagnetic excitations are 
handled by graphene. The chiral symmetry for the graphene quasiparticles in a given K-valley is either parallel 
or antiparallel to the direction of motion of the electrons and holes. An immediate consequence is that in the 
ksT Huj limit the optical conductivity for undoped graphene is independent of any materials properties, Uo = 
e^/(4h) (Kuzmenko et al. 2008). Nevertheless, even at room temperature, it is found experimentally that the optical 


absorption is very small ~ 2.3% and it depends only on the fine structure constant, a = 1/137 (Nair et al. 2008). 


Doping and gating influence the optical properties significantly leading to Pauli blocking for photons with energy 
less than 2Ef {Ep is the Fermi level) and achieving carrier concentrations which can modify the transmission in the 
visible spectrum (Li et al. 2008[ Liu et al. 2011a). Graphene plasmonics is also quite interesting since graphene 


plasmons are tunable by gating and doping and they are temperature dependent (Gangadharaiah et al. 2008). In 


addition to longitudinal plasmons, graphene can support a transverse plasmon mode. It is also interesting to note 
that the longitudinal modes are gapless, however the transverse ones exist in the window 1.7 < hw/Ep < 2 and can 
be tunable from radio to infrared frequency by doping and electric fields (Mikhailov and Ziegler, 2007). 

Plasmons in GNRs can exist in the near-infrared to far-infrared range and further tunability via a gate voltage can 


be achieved (Freitag et al. 2013). The optical properties of GNTs are also quite unique. The CNT optical activity. 


such as electron-energy-loss spectroscopy (EELS) spectra and circular dichroism, are chirality dependent (jPresselhaus 


et al. 2007 Wang et aH|2005 ). Competing effects due to Coulomb interactions and an attractive e-h coupling are 


strong in CNTs, which has led to the realization that strong excitonic effects need to be taken into account (especially 
in small diameter nanotubes) to achieve agreement with experimental optical data (Spataru et al ., 2001[ 2004 2005| ) . 

The optical response of graphene can be described by considering its 2D conductivity tensor calculated within 
the Kubo formalism (Falkovsky and Varlamov, 2007). Evaluating this general expression for the lowest conduction 
and highest valence energy bands in the q —>■ 0 approximation leads to the intraband (intra) and interband (inter) 
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contributions 




In 2 


^inter 

JA 


K^ttPuj 

2'’ sinh(;3£;) 


ln(cosh(/3A) + cosh(/i/3)) — 


2A2 


dE 


sinh(/3i?) 


Trtfiu} 7a cosh(/r/3) + cosh(^£^) 


2a;A2 


dE 


sinh(/3i?) 


1 


(11) 


cosh(/r/l) + cosh(/3i?) {Huj)^ + AE^ tt 7a cosh(/r/l) + cosh(/?i?) {Huj)^ + AE^ 
where A is an energy gap in the graphene spectrum. When A = /i = 0 and ksT <C hu!, a acquires the universal value 


when spatial dispersion is not taken into account and the difference between a^x and 

ayy (graphene is in cry-plane) is 

mostly pronounced for larger q ( 

Drosdoff et al. 

2012 

Falkovsky and Varlamov 

2007 

)■ 


The optical response properties can also be characterized by considering the longitudinal polarization function 
Xi{q,ito), which corresponds to the longitudinal component of the conductivity - (7{q,iuj) = for q —>■ 

0. Alternatively, the transverse electric (TE) and transverse magnetic (TM) excitations can be captured via the 
polarization tensor 11 calculated by a (2+1) Dirac model ( Bordag et al.\ 2009a[ Fialkovsky et al. 2011| Klimchitskaya 
et al.. 2014 Sernelius 20151. It is found that the longitudinal polarization function is related to the Boo component, 
Xi = while the transverse polarization function is Xtr = - gflloo). 

The optical response of the quasi-lD structures, such as GNRs and CNTs, follows from the Kubo formalism for 
graphene. Taking into account the zone-folded wave functions and chirality dependent energies leads to the intra and 


interband optical conductivity spectra of CNTs (Tasaki et al. 1998). Similarly, incorporating the edge dependent wave 


functions with the appropriate TB energies results in the intra and interband conductivities of zigzag and armchair 
GNRs ( ]Brey and Fertig 2006 Sasaki et al. 2011). 


C. Casimir interactions and graphene nanostructures 

1. Graphene 

The vdW/Casimir interaction involving 2D graphene can be calculated using the fully retarded expression in Eq. 
1^ with response properties (Eqs. 


11 12) corresponding to its low-energy Dirac spectrum. The boundary conditions 


are contained in the matrices Ki, 2 , whose non-zero diagonal components for a graphene/semi-infinite medium system 
reflecting the TE (ss) and TM (pp) modes are 





2'Kuiaqc^ 

1 + 2TTUjaqc^ 


. n(pp) _ 
, Ki — 


n(ss) _ 

Jto — 


fj,{iu})q - k ^ ^(pp) 
IJ,{i(jj)q + k’ ^ 


2Traq/uj 
1 + 2Traq/uj 
e(iio)q — k 
e(iuj)q + k ’ 


(13) 

(14) 


where q = + {eo/cY and k = y'/cy + p{iLo)e{iuj)uj^I. The dielectric and magnetic response functions for the 

semi-infinite medium are e{iu}) and fJ.(iu}), respectively. The reflection coefficients here are expressed in terms of the 
graphene conductivity cr, however, these can be given equivalently via other response characteristics using the relations 
discussed above. For a graphene/graphene system, the components of the K 2 matrix are replaced by the components 
of the Ml matrix. 

One of the first studies of Gasimir interactions for graphene was reported in ( [Bordag et al. 20061, where the authors 
considered graphene/perfect metallic semi-infinite medium and atom/graphene systems. The graphene is modeled as 
a plasma sheet leading to results strongly dependent on the plasma frequency. A more suitable representation of the 
graphene sheet was later considered by a series of authors by taking into account the Dirac-like nature of the carriers 
explicitly. It is obtained that the Casimir force is quite weak compared to the one for perfect metals and that it is 


strongly dependent upon the Dirac mass parameter (Bordag et al. , 2009a). Describing the graphene response via the 


2D universal graphene conductivity cto valid in the ksT ^ huj limit (Falkovsky and Varlamov 2007 Nair et al. 2008), 


other authors (Drosdoff and Woods 2009) find a unique form of the graphene/graphene Casimir force per unit area 

'-o"T-' ' -' 

This result shows that the distance dependence is the same as the one for perfect metals 
= — 24 odi 1 however the magnitude is much reduced, Eq/E„ 


A Eo — 3 h< 7 o _ _ 

A 87Td* 32'n-d-* 


whose Casimir force is ^ 


0.00538. It is interesting to 

note that retardation does not affect the interaction (no speed of light c) and k is canceled after taking into account 
that CTo = Y!Ah. 
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The graphene interaction has also been investigated via the non-retarded Lifshitz formalism in Eq. ([Dobson et al. 


2006 Gomez-Santos 

2009 

Sarabadani et al. 

2011 

Sernelius 

2015 


where the polarization and Coulomb interaction 
are calculated with the RPA approach. The RPA is a useful tool to study long-ranged dispersive interactions as it gives 


a natural way to take into account the electron correlation effects of each object and spatial dispersion (Dobson, 2011 


Fetter and Walecka 1971), as discussed earlier. It has been found that for separations d > 50 nm the non-retarded 
Lifshitz approach results in a graphene/graphene force of the form ^ — Jl, where the magnitude of the constant B 


agrees with the results from the retarded Casimir calculations (Drosdoff et al. 2012 Drosdoff and Woods 2009). It 


is thus concluded that the graphene/graphene interaction is determined by the non-retarded TM mode contribution 
(captured in the longitudinal polarization) even at distances corresponding to the Casimir regime. These results are 
truly remarkable since the vdW/Casimir interaction appears to be independent of all of the graphene properties in the 
low T and/or d > 50 nm regime. A further interpretation can be given by noting that the electromagnetic fluctuations 
exchange occurs at speed vp (Eq. 10) rather than the speed of light. This means that the typical thermal wavelength 
At = hc/ksT, which sets the scale where quantum mechanical (d < At) or thermal (d > At) fluctuations dominate 
the interaction, becomes At = hvp/kpT. The quantum mechanical contributions determine the graphene interaction 
at separations d < At ~ 50 nm as opposed to d < At 7 /rm for typical metals and dielectrics at T ~ 300 K. 
The thermal fluc tuations for graphene become relevant at much reduced distances, and for d > At the interaction is 


Ft 


C(3) fcflT 

Stt 


(Gomez-Santos 


2009) where /'(n) is the Riemann zeta-function. Essentially, vp takes the role of the 


speed of light enhancing the importance of the zero Matsubara frequency at much lower T and smaller d as compared 
to conventional metals and dielectrics. 

For closer separations (d < 50 nm), a more complete model for the graphene properties is needed. Deviations from 
the asymptotic behavior at low T are found (Drosdoff and Woods 2009) by using the graphene optical conductivity 
taken into account by a Drude-Lorentz model that corresponds to higher frequency range tt —> tt* and cr —> cr* 
transitions. Recently, it has been shown that the Casimir interaction in a stack of identical grap hene layers exhibits a 
fractional distance dependence in the energy {E ~ d“®/^) as a result of the Lorentz oscillators (Khusnutdinov et al. 


2015). Other researchers ( Gould et al.\ 2013a[ 2009, 2013b[[2008 Lebegue et al. 2010) have utilized the RPA approach 


combined with first-principles calculations for the electronic structure to investigate the non-retarded interaction at 
very short separations (d < 10 nm) for infinite number of parallel graphene layers. Interestingly, the interaction 
energy is found to be E ^ d“^. This insulator-like behavior is attributed to the full energy band structure (beyond 
the two-band model in Eq. 10) and the associated higher transitions in the response properties. Other authors 


(Sarabadani et al. 2011) have also considered the vdW interaction in a multi-layered graphene configuration within 


the RPA, however, the reported unusual asymptotic distance dependences may be an artifact of the considered finite 
graphene thickness. 

It has also been shown that temperature together with other factors, such as doping or external fields, affect the 
graphene thermal and quantum mechanical regimes in an intricate way. In particular, the classical Casimir/vdW 
interaction determined by the thermal fluctuations has been examined by several authors in different situations. 


(Fialkovsky et al. 


2011) have used the polarization tensor and corresponding reflection coefficients to e xpress the 
dominating thermal fluctuations regime in terms of the fine structure constant a as La™ ^ 


2C(3) 


(Sernelius 


2011) has utilized the longitudinal graphene response in Eq. [^to show that doping plays an important role in the 


interaction at larger separations, as the force can be increased by an order of magnitude for large degrees of doping. 
( Bordag et al.\ 2012| Klimchitskaya and Mostepanenko 2013| ) have used the fully relativistic Dirac model with the 
T-dependent polarization tensor to investigate how a finite mass gap A in the Dirac model affects these regimes in 
graphene/graphene and graphene/dielectrics. It is found that for kpT A, thermal fluctuations are not important, 
while for A < kpT the thermal effects become significant, as shown on Fig. [^a, b). The thermal and quantum 
mechanical regimes have also been studied by (Drosdoff et al. 2012) via the longitudinal thermal conductivity, which 
includes spatial dispersion, an energy gap and chemical potential in the Dirac model. These authors have shown that 
tuning A and /i can be effective ways to modulate the interaction, however, the spatial dispersion does not play a 
significant role except for the case of small A and low T, as shown in Fig. id, e). 

Recent studies (Bordag et al. 2015 Klimchitskaya and Mostepanenko] |2014[ [Klimchitskaya et ^ 2014 Sernelius 


2015) provide a thorough analysis of the balance between the thermal and quantum mechanical effects in the graphene 


Gasimir interaction. It is shown that equivalent representations within the temperature dependent longitudinal and 
transverse polarization and the temperature dependent polarization tensor are possible. The comparison between 
results from the temperature dependent polarization tensor and density-density correlation function show that at low 
T, both approaches give practically the same results proving that retardation and TE polarization are unimportant. 
For T 7 ^ 0, deviations are found as shown in Fig. [^c) for graphene/graphene and graphene/metal configurations. 


A paper by (Dobson et al. 2014) reveals that the collective excitations beyond the RPA approximation may be 
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- A = 0.1 eV 



A (eV) n (eV) 


FIG. 5 (Color Online) (a) Casimir pressure between two graphene sheets at separation d = 30 nm as a function of temperature 
for different values of the gap A (Figure adapted from (Klimchitskaya and Mostepanenko 20131); (b) Relative thermal correction 
SF{%) — {F{T) — F{T = 0))/F{T = 0) for the graphene-Si plate interaction at separation d = 100 nm (Figure adapted from 
( Bordag et al.\ 20121). (c) Relative deviation 5F{%) = {Fdd{T) — Fpt{T))/Fpt{T) for graphene-graphene (red) and graphene-Au 
plate (blue) interactions, where Fdd{T) is the Casimir force calculated via the density-density correlation function and Fpt{T) is 


the Casimir force calculated via the polarization tensor (Figure adapted from (Klimchitskaya et al. 
graphene force normalized to Fo = —3e^/(327rd^) with and without spatial dispersion in the grap 


as a function of: (d) the gap A and (e) the chemical potential /r (Figures from ([Drosdoff et al.\ 


20141). Casimir graphene- 


lene conductivity at T = 0 

2 ^). 


quite important, qualitatively and quantitatively, for the graphene non-retarded vdW interaction. In general, it is 
assumed that higher vertex corrections may change the magnitude of the force somewhat, but not the asymptotic 
distance dependence. However, this may not be the case for graphene as the type of renormalization yields very 
different results. The renormalization-group method (Kotov et al. 2012 Sodemann and Fogler 20121 results in a 


weak correction to the interaction energy as opposed to the two-loop level in the large-iV limit approach (Das Sarma 


et al. 

2007), where the characteristic distance dependence has a different power law ( 

Dobson et al. 

2014 


These 

findings indicate that graphene may be the first type of material for which RPA is not enough to capture the vdW 
force in the quantum limit. Along the same lines, (Sharma et al. 2014) have shown that for strained graphenes, where 
electron-electron correlations beyond RPA are much more pronounced, corrections to the vdW interaction, consistent 
with the renormalization-group model, are found. 

Besides the fundamental questions regarding basic properties of graphene Casimir/vdW interactions, other and more 
exotic applications of this phenomenon have been proposed. For example, (Phan et al ., |2012' ) have proposed that a 
graphene flake suspended in a fluid, such as Teflon or bromobenzene, can serve as a prototype system for measuring 
thermal effects in Casimir interactions. The balance of gravity, buoyancy, and the Casimir force on the flake creates 
a harmonic-like potential, which causes the flake to be trapped. By measuring changes in the temperature dependent 
frequency of oscillations, one can potentially relate these changes to the Casimir interaction. Alternative ways to tailor 
the graphene Casimir interaction have also been recognized. For example, (Svetovoy et al. 2011) have shown that 
the thermal effects can be enhanced or inhibited if one considers the force between graphene and different substrates. 


(Sernelius 2012) finds that retardation due to the finite speed of light can also be made prominent depending on 


the type of substrates graphene interacts with. (Drosdoff and Woods, 2011) have proposed that metamaterials with 
magnetically active components can result in a repulsive Casimir force. (Phan et al. 2013) have shown a regime where 
repulsion can be achieved with a lipid membrane. Also, Dirac carriers with constant optical conductivity result in 
unusual Casimir effects behavior in nonplanar objects. For example, the interaction on a spherical shell with a =const 
has markedly different asymptotic behavior and sign when compared to the one for a plasma shell or for planar sheets 


with a =const (Bordag and Khusnutdinov 2008 Khusnutdinov et al. 2014). 


The Casimir-Polder force involving atoms and graphene sheets has also been of interest. Such studies are relevant 
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not only fundamentally, but also for other phenomena, including trapping or coherently manipulating ultra-cold atoms 
by laser light (Bajcsy et al.. 2009 Goban et al.\ 2012 Ito et al. 1996). The theoretical description follows from Eg. 


by considering one of the substrates as a rarefied dielectric (Dzyaloshinskii et al. 1961 Lifshitz and Pitaevskii 1980 


Milonni 1993). In addition to atom/graphene (Judd et al.. 20111, configurations containing additional substrates have 
been studied ( Chaichian et aL[|2012 ). Authors have suggested that it may be possible to observe quantum reflection 
of He and Na atoms via Casimir-Polder interaction as means to discriminate between the Dirac and hydrodynamic 


model description for graphene (Churkin et al. 2012). Casimir-Polder shifts of anisotropic atoms near multi-layered 


graphene sheets in the presence of a Huttner-Barnett dielectric (a linearly polarizable medium, which is modeled 
by microscopic harmonic fields), have also been calculated (Eberlein and Zietal, 2012). Thermal fluctuation effects 
in atom/graphene configurations can also be much stronger due to the reduced thermal wavelength Ay. Thermal 


the zero Matsubara frequency giving rise to Ft = — 


et al. 


et al. 2012 Kaur et al. 2014). Interesting possibilities for temporal changes in the atomic spectrum affecting the 


2015 

Chaichian et al. 

2012 

Drosdoff 


graphene sheet by creating ripples have also been suggested (Ribeiro and Scheel 2013a). The Casimir-Polder potential 


has further been explored possibilities for shielding vacuum fluctuations in the framework of the Dirac model (Ribeiro 


and Scheel 2013b). 


2. Quasi-ID graphene nanostructures 


Investigating atom/CNT interactions is of utmost importance for applications, such as trapping cold atoms near 
surfaces (Goodsell et al.\ 2010 Petrov et al. 2009), manipulating atoms near surfaces for quantum information 


processing (Schmiedmayer et al.\ 2002), and hydrogen storage (Dillon et al. 1997). CNT/CNT interactions are 


relevant for the stability and growth processes of nanotube composites (Charlier et al.\ 2007). To calculate the 


interaction, one must take into account the cylindrical boundary conditions. Researchers have utilized scattering 
techniques to study the distance dependence involving metallic wires with Dirichlet, Neumann, and perfect metal 


boundary conditions ( 

Emig et al. 

2006 

Noruzifar et al. 

2011). Inclined metallic wires have also been considered 

(Dobson et al. 

2009 Noruzifar et al. 20E 

). Calculations for CNT interactions, however, are challenging as one has to 


take into account simultaneously the chirality dependent response properties and the cylindrical boundary conditions 
for the electromagnetic fields. 

The Lifshitz approach has been applied to CNTs via the proximity force approximation, which is typically appropri¬ 
ate at sufficiently close separations (Blocki et al ., 1977| ) . The cylindrical surface is represented by an infinite number 
of plane strips of infinitesimal width, which are then summed up to recover the CNT surface. This method has been 
applied to atom/single-walled nanotube and atom/multiwall nanotube treated as a cylindrical shell of finite thickness 
(Blagov et al. 2005, 2007| Bordag et al. 2006[ Churkin et al.. 2011 Klimchitskaya et al. 2008). In these studies, 
the dielectric response of the nanotubes is not chirality dependent. (Blagov et al. 2005) uses extrapolated dielectric 


Function for graphite. 

(Bordag et al. 

2006 

) uses the response to be due to a surface density of the 7r-electrons smeared 

over the surface. 

(Blagov et a/.[|2007h utilizes a free-electron gas representation for the cylindrical CNT surfaces, while 

(Churkin et al. 

2011) 

takes the graphene Dirac and hydrodynamic models. In these works, the interaction energy is 


always of the form E = —C 3 {d)/(P, where the coefficient Cs^d) is also dependent on the cylindrical curvature and 
atomic polarizability. 

Interactions between nanotubes in a double-wall configuration have been calculated using the QED approach suitable 
for dispersive and absorbing media as well (Buhmann and Welsch, 2007). Within this formalism the boundary 


conditions are taken into account by solving the Fourier domain operator Maxwell equations using a dyadic Green’s 
function, which also allows the inclusion of the chirality dependent response properties of the individual nanotubes. 
The calculations utilize the fluctuation-dissipation theorem and the force per unit area is the electromagnetic pressure 
on each surface expressed in terms of the Maxwell stress tensor ( [Tomas 2002). For planar systems, the QED and 
the Casimir/Lifshitz theory lead to the same expression (Eq. [^, which has also been shown for systems involving 
graphene (Drosdoff and Woods, 2009 Hanson 2008). The QED method, applied to the interaction in various double- 


walled CNTs, revealed that the chirality dependent low-energy surface-plasmon excitations play a decisive role in 


the interaction (Popescu et al. 2011 Woods et al. 2011). The attractive force is actually dominated by low-energy 


inter-band plasmon excitations of both nanotubes. The key feature for the strongest attraction is for the CNTs to 
have overlapping strong plasmon peaks in the electron energy localization spectra (EELS). This is true for concentric 
(n, n) armchair CNTs, which exhibit the strongest interaction as compared to tubes with comparable radii, but having 
other chiralities, as shown in Fig. [^ The results are consistent with electron diffraction measurements showing that 
the most probable double-walled CNT is the one in which both tubes are of armchair type ([Harihara et al. 2006). 


































































































































































































18 



FIG. 6 (Color Online) Electromagnetic pressure on each nanotube in a double-wall CNT system as a function of separation. 
The inset shows calculated EELS for several armchair (n, n) and zigzag (n, 0) nanotubes as a function of frequency in eV. The 
attraction is strongest between two concentric armchair nanotubes due to the presence of strong overlapping low frequency 
peaks in the spectra. The notation (m, 0)@(n, n) corresponds to (m, 0) as the inner tube and (n,n) as the outer tube. 


This indicates that the mutual Casimir interaction influenced by the collective excitations may be a potential reason 
for this preferential formation. 

The Casimir-Polder interaction involving CNTs has also been considered via the QED formalism. For this purpose, 
one utilizes a generalized atomic polarizability tensor containing dipolar and multipolar contributions and a scattering 


Green’s function tensor expressed in cylindrical wave functions (Buhmann et aZ. 2004 Li et al. 2000 Tai 1994). It has 


been shown that the chirality dependent CNT dielectric function plays an important role determining the strength of 


the atom/nanotube coupling (Fernari et al. 2007). The QFD approach has also been used initially by (Bondarev and 
Lambin[ 2004[ 20051, where the nonretarded interaction potential is equivalently given in terms of photonic density 


of states. These studies also show that the interaction is sensitive to the CNT chirality (Rajter et al. 2013, 2007). It 


was found that the stronger optical absorption by the metallic CNTs suppresses the vdW atomic attraction, which 
can be of importance to tailor atomic spontaneous decay near CNT surfaces. 

For GNRs the situation is even more technically difficult as compared to nanotubes since analytical results for 
the boundary conditions for strip-like systems are not available. Nevertheless, a perturbation series expansion of the 
Lifshitz formula in the dilute limit separates the geometrical and dielectric response properties contributions into 
convenient factor terms, which can be quite useful to study the finite extensions of the nanoribbons in their vdW 


interaction (Stedman et al. 2014). A recent study has also shown that a non-retarded Lifshitz-like formula for the 


vdW interaction between parallel quasi-ID systems having width W <s:. d can also be derived (Drosdoff and Woods 


2014). The force per unit length is written in terms of a TM-like ’’reflection coefficient” containing the CNR response 


properties, which makes the expression reminiscent of the Lifshitz vdW expression for planar objects in Fq. This 
is quite appealing as it presents a general way to calculate vdW interactions in any type of ID parallel systems. 


2014). 


Applying this theory to GNRs described by their specific response properties ( ]Brey and Fertig 2006) shows that 
the chemical potential is crucial in the interplay between quantum mechanical and thermal effects in the interaction. 
A /i-dependent transition between these two regimes is reported correlating with the onset of intraband transitions 

While GNRs with p, = 0 behave like typical dielectric materials with a vdW force 
g {Eg is the energy gap in the CNR band structure), the interaction becomes completely 
thermal with a characteristic behavior F = ~ • For semiconductors, such as GaAs wires however, the 

thermal fluctuations dominate the interaction completely. This is at complete odds with the dispersive interaction 
involving 2D or 3D materials, where thermal effects are typically very small (Klimchitskaya et al. 2009). It turns 


(Drosdoff and Woods 
F ~ when [j, 


E„ 


out that for GaAs quantum wires the plasma frequency is much reduced as compared to 3D GaAs (Das Sarma et al. 


1996) resulting in the force being dominated by the n = 0 Matsubara term. 




























































19 


D. Materials with topologically nontrivial phases 


In addition to graphene, there are other materials with Dirac spectra and topological insulators (TIs) have a special 
place in this class of systems. TIs are a new phase of matter with nontrivial topological invariants in the bulk electronic 
wave function space. The topological invariants are quantities that do not change under continuous deformation, and 


they lead to a bulk insulating behavior and gapless surface Dirac states in the band structure (Ando 2013 Cayssol 


2013 

Hasan and Kane 

o 

o 

Qi et al. 2008 

Qi and Zhang 

2011 


realization that a strong spin-orbit coupling can result in a TI phase with several materials being proposed as possible 
candidates, including BixSbi_x, Bi 2 Se 3 , Bi 2 Te 3 , TlBiSe 2 among others (Ando 2013 Chen et al. 20091. In the low 


momentum limit the 2D states, which are topologically protected by symmetry invariance, are described by a helical 


version of the massless Dirac Hamiltonian (Liu et al. 2010 Zhang et al. 2009) 


Hsurf — Z • ((J X k) , 


(15) 


where z is the unit vector perpendicular to the surface (located in the a;y-plane), a = (cTa,, try, Uz) are the Pauli 
matrices and k is the 3D wave vector. 


Topologically nontrivial materials can be classified via their symmetries and dimensions (Hasan and Kane, 2010 


Kitaev, 2009 Ryu and Takayanagi 2010 Schnyder et al. 2008) or by dimensional reduction (Qi et al. 2008). Three¬ 


dimensional systems are characterized by time reversal (TR) symmetry leading to each eigenstate of the above 
Hamiltonian being accompanied by its TR conjugate or Kramers partner (Schnyder et al.. 2008). Experimentally, 
however, one observes an odd number of Dirac states. This is understood by realizing that the Dirac cones of the 
Kramer’s pairs appear on each side of the surface of the material and the cone in empty space cannot be detected. In 
addition, these surface states are protected from backscattering by the TR symmetry, which makes them insensitive 
to spin-independent scattering - a useful feature for quantum computation applications (Leek et al. 2007). Chern 


Insulators (CIs) are essentially two-dimensional TIs and their low-energy band structure, also described by Eq. (15), 


consists of even number of helical edge states. CIs have strong enough interband exchange energy, responsible for the 
hybridization of the surface states from the Dirac cone doublets. CIs states are further described by a topological 
integer Chern number C S Z quantified as C = ^ sign(Ai), where N denotes the (even) number of Dirac cones 
and Ai is the mass gap of each Dirac cone. The mass gap can be tailored by an applied magnetic field or other means 
and it can be positive or negative. 

The properties affecting the vdW/Casimir interactions in systems with topologically non-trivial phases are linked to 
the dimensionality and response characteristics of the involved Dirac materials. The optical conductivity components 
of TIs and CIs involving the low-energy Dirac carriers have been obtained within the standard Kubo approach or 
the quantum kinetics equation method (Rodriguez-Lopez and Grushin 2014 Tse and MacDonald 2010| . Analytical 
representations for the longitudinal surface optical conductivity in the small temperature regime ksT min(|/rF|, | A|) 
(here fip is the Fermi energy relative to the Dirac point) and small disorder have been found (Chen and Wan 2011[ 
2012[ Grushin et al.. 2012 Tse and MacDonald[ 2010[ 2011) with expressions similar to the ones for graphene (Eqs. 
TTf 12). In addition, the topologically protected surface states lead to a strong quantum Hall effect (QHE) without an 
external magnetic field whenever perturbations breaking the TR symmetry induce a gap in the band structure. For 


the low-energy carriers in Eq. 15 the associated surface Hall conductivity has the following expression at imaginary 
frequency 


(7xy{iuj) = - 


acA 

2'Khuj 


tan 


-1 


Huj \ 1 / huj 

— - tan"^ 

2eJ V2|A| 


(16) 


Here Cc is the energy cutoff of the Dirac Hamiltonian, which we associate with the separation between the Dirac point 
and the closest bulk band. 

For the Gasimir interaction involving GIs, it is important to note that these materials can exhibit a quantum 
anomalous Hall effect at zero frequency or in the absence of an external magnetic field. By tuning the mass gap 
(via doping or changing the magnetization of the involved material), one can eliminate the 2D optical conductivity 
cTxxi^ = 0, \hf\ < |A|) = 0, while the Hall conductivity becomes axy{uj = 0, |^_f| < |A|) = ^sign(A). After summing 
up the contributions from all Dirac cones, one obtains a quantized Hall conductivity in terms of the Ghern number 


^xy 


iu = 0) = if C. 


Inducing a mass gap has important consequences for the surface Hall response in 3D TIs, as well. By applying 
an external magnetic field, it is possible to realize the fractional quantum Hall effe c t with a quantized conductivity 
axy{oJ = 0, |A| > Ifip]) = ^ (I +^)j where n is an integer (Hasan and Kane 2010| Zheng and Ando 2002). 
Nevertheless, one also needs to add the bulk dielectric response! Typically, a standard Drude-Lorentz model is 
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FIG. 7 (Color Online) (a) Energy band structure of a Cl calculated via a generic two-band tight binding model ( Grushin et al.\ 

(b) Real 
correspond 


20121. The low-momentum limit of the energy band structure is consistent with the Dirac Hamiltonian in Eg. |15| 
part of axy{u) (left panel) and a^yiiuj) (right panel); (c) Real part of (Jxx{oo) (left panel) and axx(iu>) (right panel) cc 
to the energy band structure from (a). The calculations are performed with C = 1, A = 0.25t and e = 2.25t {t denotes the 
hopping integral for the employed lattice model; here it is taken to be equal to the frequency bandwidth). The conductivities 
are in units of ac/( 27 r). The analytically found and a^yj^) are in excellent agreement with numerically evaluated Kubo 


expressions. (Figure taken from (Rodriguez-Lopez and Grushin 20141) 


sufficient, and authors have shown that specifically for the Casimir interaction a single oscillator for the dielectric 


2011). Therefore, the bulk response can be considered as e{iuj) = eo + 
or and ujr is the location of the resonance. 




where We is the strength of the 


function is enough to capture the characteristic behavior (Chen and Wan, 2011 Grushin and de Juan 2012 Grushin 

-^^ 

et al. 
oscilla 

The surface response properties dramatically affect the electrodynamics in topologically nontrivial materials in 3D. 
In fact, the electrodynamic interaction can be described via generalized Maxwell equations containing a magnetoelec¬ 
tric coupling due to the surface Hall conductivity. Equivalently, this generalized electrodynamics includes an axion 
field 0{r,t) manifested in a Chern-Simmons term in the Lagrangian, Cg = E • B, whose role is to preserve the 

TR symmetry in the Maxwell equations (Wilczek, 1987). While 0{r,t) depends on position and time in general, for 
topological insulators, this is a constant field, such that d 7 ^ 0 in the bulk and 0 = 0 in the vacuum above the surface 
of the material. We further note that the quantization of the Hall effect in 3D TIs is inherited in the axion term 
according to 9 = {2n + l) 7 r 


Essin et al. 

2009 

Qi et al. 

2008 


The concept of axion electrodynamics was first proposed in high energy physics as a possible means to explain dark 
matter (Peccei and Quinn 1977 fWilczek 1987), and now axion-type of electromagnetic interactions appear in the 
description of condensed matter materials, such as TIs. The axion field originating from the topologically nontrivial 
surface states leads to many new properties, including induced magnetic monopoles, quantized Faraday angle in 
multiple integers of the fine structure constant, and large Kerr angle (Qi et al. 2009 Tse and MacDonald, 2010 


2011 Wilczek, 1987). The modified electrodynamics due to the Ghern-Simmons term with the associated boundary 


conditions is also of importance to the Gasimir interaction, as shown earlier from a high energy physics perspective 
(Bordag and Vassilevich, 2000). 


E. Possibility of Casimir repulsion in topological materials 


The underlying electronic structure of the materials and their unconventional Hall response, however, open up 
opportunities to explore the Gasimir effect in new directions. Fig. [^a) depicts the low-energy Dirac band structure 
for an appropriate lattice model for a GI with the associated longitudinal (axx) and Hall (axy) conductivities (Grushin 
et al. 2012 Rodriguez-Lopez and Grushin 2014). The reflection matrices in the Lifshitz expression from Eq. |^have 
been determined for two semi-infinite TI substrates with isotropic surface conductivity aij and dielectric and magnetic 
bulk response properties taken as diagonal 3D matrices e = e(tLi)I and /r = respectively. Generalizations due to 

non-local effects and anisotropies in the response (Grushin et ail 2011 1 , as well as finite width substrates (Peterson 


and Strom 1974), can also be included. The reflection coefficients for GIs follow from the ones for the 3D TIs simply 
by setting e,/i — )■ 1 (Martinez and Jalil, 2013 Rodriguez-Lopez and Grushin 2014 Tse and MacDonald 2012). It 
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FIG. 8 (Color Online) (a) Casimir interaction energy E{d) between two CIs in u nits of Eqjd) = —hca^/(8n^dr) as a fu nction 


A phase diagram for the interaction energy between two semi-infinite TIs substrates for all separation scales. 


2014| )). (b) 
The 


of d = td/{hc) for different values of Ci 2 and |A| = t. (Figure adapted from (Rodriguez-Lopez and Grushin 

’ - 

2n > 

parameters cue and cjh correspond to the strength and location of the Drude-Lorentz oscillator, respectively. The repulsion for 
large separations at T = 0 if is given by the red region (outlined by the red line), while the repulsion for all separations at high 
T is given by the green region (outlined by the green line). Here, ^ = (’1-+ |) S'lid Q\ = —02 = d for the substrates. Repulsion 
is observed for a large range of Chern numbers. 


turns out, however, that in all cases the surface Hall conductivity is a key component in understanding the asymptotic 
distance dependence, magnitude, and sign of the interaction. 

Fig. j^a) summarizes results for calculated Casimir energies at the quantum mechanical regime (low T and/or large 
d). The graph indicates that there is a change of distance dependence behavior when comparing the small and large 
d asymptotics for interacting CIs. The analytical expressions for the conductivity components, which agree very well 
with the numerical Kubo formalism calculations according to Fig. Hb, c), are especially useful in better understanding 
of the underlying physics of the Casimir energy. It has been obtained that the energy at small d is determined by the 
longitudinal component of the conductivity and the interaction is always attractive |Rodriguez-Lopez and Grushin 


20141. For large d, however, it is possible to achieve repulsion if the two CIs have Chern numbers C 1 C 2 <0 and the 


Hall conductivity is much larger than the longitudinal one. The interaction energy in this case is a non-monotonic 
function of distance and it is quantized according to E ^ CiC 2 , which also indicates that the strongest repulsion 
occurs for materials with large Chern numbers. It is concluded that if repulsion is desired, one must search for Cl 
materials with vanishing diagonal conductivity components and strong Hall conductivity capable of sustaining much 
enhanced C numbers. Let us note that a quantized Casimir interaction may be typical for materials that can support 
a strong Hall effect. In fact, such phenomenon was predicted to exist in graphene/graphene systems at sufficiently 
large separations with an external magnetic field (Tse and MacDonald, 2012). The associated Landau-level filling 


factors in the Hall conductivity lead to a quantization condition in the force, and similar findings have been reported 
for atom/graphene configurations (Cysne et al.. 2014). 

The Casimir interaction between TIs has also been studied, in which case the bulk dielectric response is included 
(Sec. HI.D). Recent work (Grushin and Cortijo 2011) has shown that the energy has unique characteristics due to 
the balance between the bulk and surface states contributions mediated by the axion term 6 . Fig. ib) summa¬ 
rizes numerical results for the interaction energy phase diagram showing repulsive and attractive regimes depending 
on the Drude-Lorentz parameters (^) and the surface contribution ( 6 ). Reported analytical calculations enable a 
better understanding of the important factors determining the interaction in various limits. It is found that if the 
bulk response is treated as a single Lorentz oscillator, it is possible to obtain Casimir repulsion at larger separa¬ 
tions (Rodriguez-Lopez, 2011), where the surface contributions through the Hall conductivity dominate the response 
provided 0i02 < 0. For shorter separations the bulk response becomes dominant, and the Casimir interaction is 
attractive. One also obtains attraction at all distance scales when di02 > 0 (Rodriguez-Lopez, 2011). Other authors 


have predicted that repulsion is also possible in the regime of short separations, however, this is considered to be an 


artifact of a frequency independent surface conductivity taken in the calculations (Chen and Wan 2011 Grushin and 


Cortijo 

2011 

Grushin et al. 

2011 

Nie et al. 

2013 


2013). In addition, several recent works have shown that the behavior 


of the Casimir interaction and the existence of repulsion, in particular, depend strongly on the magnitude of the finite 


mass gap, the applied external magnetic field, and the thickness of the TI slabs ( 

Chen and War 

2011 

Nie et al. 

2013). 

The thermal Casimir interaction between TIs has also been studied (Grushin and Cortijo 

2011 

Rodriguez-Lopez 


2011). The energy corresponding to the n = 0 Matsubara term depends strongly on the axion fields. The interaction 


is found to be attractive when 6162 > 0. However, thermal Casimir repulsion is obtained at all distances for 6 i 02 < 0 
as shown in Fig. |^b). Similar considerations for repulsion may apply for CIs, since their thermal Casimir energy can 
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be obtained analogously to be proportional to the surface Hall conductivities . These results indicate that 

topological Dirac materials, such as 2D CIs and 3D TIs, may be good candidates to search for a repulsive thermal 
Casimir interaction. 

The long-ranged dispersive interactions involving systems with nontrivial topological texture are complex phenom¬ 
ena. Materials with Dirac carriers lend themselves as templates where concepts, typically utilized in high energy 
physics, cross over to condensed matter physics with vdW/Casimir interactions as a connecting link. Topologically 
nontrivial features in the electronic structure and optical response properties result in unusual asymptotic distance 
dependences, an enhanced role of thermal fluctuations at all distance scales, and new possibilities of Casimir repul¬ 
sion. Ongoing work in the area of Dirac materials will certainly continue stimulating further progress in the field of 
vdW/Casimir physics and further widening the scope of fluctuation-induced phenomena. 


IV. STRUCTURED MATERIALS 


The geometry of the interacting objects and the interplay with the properties of the materials is also of interest for 
tailoring the Casimir force. Structured materials, including metamaterials, photonic crystals, and plasmonic nanos¬ 
tructures, allow the engineering of the optical density of states by proper design of their individual components. As a 
result, one is able to manipulate the interaction utilizing complex, non-planar geometries. Recent experimental studies 


have begun the exploration of such geometry effects particularly with dielectric and metallic gratings (Chan et al. 


2008 Intravaia et al. 2013). The theoretical description has been challenging due to the complex dependence of dis¬ 


persive interactions upon non-planar boundary conditions. One approach relies on effective medium approximations, 
where the emphasis is on models of the dielectric response of the composite medium as a whole. The second approach 
deals with particular boundary conditions via computational techniques. While the interaction of electromagnetic 
waves with metallic and dielectric structures of complex shapes is well established in classical photonics, the main 
challenge stems from the inherently broad band nature of Casimir interactions, where fluctuations at all frequencies 
and wave-vectors have to be taken into account simultaneously. 


A. Metamaterials 


Electromagnetic metamaterials are composites consisting of conductors, semiconductors, and insulators, that res¬ 
onantly interact with light at designed frequencies. The individual components make up an ordered array with unit 
cell size much smaller than the wavelength of radiation. As a result, an electromagnetic wave impinging on the 
material responds to the overall combination of these individual scatterers as if it were an effectively homogeneous 
system. Metamaterials were speculated almost 50 years ago by Victor Veselago ( Veselagol 1968| ), who was the first 
to explore materials with negative magnetic permeability in optical ranges. However, it was over twenty years ago 
that John Pendry proposed the workhorse metamaterials’ structure, the split-ring resonator (SRR), that allowed an 
artificial magnetic response and was a key theoretical step in creating a negative index of refraction ( jPendry et al. 
1999). David Smith and colleagues were the first to experimentally demonstrate composite metamaterials, using a 


combination of plasmonic-type metal wires and an SRR array to create a negative effective permittivity eeff(<^) and 
a negative effective permeability /ieff(w) in the microwave regime (Shelby et al.. 2001). Many exotic phenomena have 
been discovered afterwards, including negative index of refraction, reversal of Snell’s law, perfect focusing with a flat 
lens, reversal of the Doppler effect and Cherenkov radiation, electromagnetic cloaking, and transformation optics. 

Such materials are of great interest for Casimir force modifications. Casimir repulsion was predicted by Boyer 


(Boyer 1974) between a perfectly conducting plate and a perfectly permeable one, but it may also occur between 


real plates as long as one is mainly (or purely) nonmagnetic and the other mainly (or purely) magnetic (Kenneth 
et al. 2002). The latter possibility has been considered unphysical (lannuzzi and Capasso, 2003), since naturally 


occurring materials do not show strong magnetic response at near-infrared or optical frequencies, corresponding to 
gaps d = 0.1 — 10/rm. However, recent progress in nanofabrication has resulted in metamaterials with magnetic 
response in the visible range of the spectrum (Shalaev 2007), fueling the hope for “quantum levitation”. 


The Casimir force for structured materials with unit cells much smaller than the wavelength of light can be calculated 
via Eq. for magnetodielectric media with reflection coefficients for two identical substrates (Mi = M 2 = M) found 
as: 


Jliss) ^ 


lies{iuj)q - k 


, R^PP^ = 


eesiiuj)q - k 
eesiioJn)q- k’ 


(17) 
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FIG. 9 (Color Online) Casimir force per unit area A between a metallic semi-space and an anisotropic metallic-based planar 
magnetic metamaterial with a weak Drude background. The filling factors f\\ and /j_, parallel and orthogonal to the vacuum- 
metamaterial interface, account for the fraction of metallic structure contained in the metamaterial. The SRR Drude parameters 
are Qd = ‘I'KcjA = 1.37 x 10^® rad/sec and 70 = 0.006 Do (corresponding to silver) and its Drude-Lorentz parameters are 
De/Do = 0.04, Dm/Do = 0.1, tJe/Do = Wm/Do = 0.1, and 7 e/Do = 7 m/D_D = 0.005. The Drude parameters for the metallic 
semi-space are D = 0.96 Do and 7 = 0.004 Do. Temperature is set to zero, and the distance between the bodies is fixed at 
d — A. Figure taken from (Rosa et al.\ 


2008a I. 


where q = + (w/c)^ and k = -|- /reffCeffW^/c^. Calculations based on this approach have suggested that 

left-handed metamaterials might lead to repulsion (Henkel and Joulain 2005 Leonhardt and Philbin 20071. Meta¬ 


materials, however, typically have narrow-band magnetic response and are anisotropic. Thus questions naturally arise 
concerning the validity of such predictions for real systems. Given that the Lifshitz formula is dominated by low- 
frequency modes w < c/d, a repulsive force is in principle possible for a passive left-handed medium as long /Xeff(fa;) 
is sufficiently larger than in that regime. Then the repulsion is a consequence of the low-frequency response 

behavior and not of the fact that the medium happens to be left-handed in a narrow band about some real resonant 
frequency. Application of the Lifshitz formalism requires the knowledge of Cegiito) and /reff(ia;) for a large range, up 
to the order of w = c/d. Such functions can be evaluated via the Kramers-Kronig relations in terms of eeff(a;) and 
/ieff(a;) at real frequencies. The point about the broad band nature of the response properties is very important, 
as it shows that knowledge of a metallic-based metamaterial near a resonance is not sufficient for the computation 
of Casimir forces: the main contribution to ees^iuj) and ^eff(*w) typically comes from frequencies lower than the 
resonance frequency. This also implies that repulsive forces, if any, are in principle possible without the requirement 
that the metamaterial resonance should be near the frequency scale defined by the inverse of the gap of the Casimir 
cavity. 

In typical metamaterial structures, the effective electric permittivity eeff(w) and magnetic permeability ^eff(w) close 
to the metamaterial resonance are well described in terms of a Drude-Lorentz model, 


eeff(w),Aieff(w) = 1 - 








(18) 


in which De(Dm) is the electric (magnetic) oscillator strength, uJeitOm) is the metamaterial electric (magnetic) res¬ 
onance frequency, and 7e(7m) is a dissipation parameter. These parameters depend mainly on the sub-wavelength 
geometry of the unit cell, which can be modeled as a LRC circuit. For metamaterials that are partially metallic, such 
as SRRs (operating in the GHz-THz range) and fishnet arrays (operating in the near-infrared or optical) away from 

resonance, it is reasonable to assume that the dielectric function also has a Drude background cdIoj) = 1 — 

(here Dd is the plasma frequency and 7 d is the Drude dissipation rate). As the Drude background clearly overwhelms 
the resonant contribution at low frequencies, it contributes substantially to the Casimir force between metallic meta¬ 
material structures. Effects of anisotropy, typical in the optical response of 3D metamaterials and of 2D metasurfaces, 
can also be incorporated (Rosa et al. 2008a|. Fig. which depicts the Casimir force between two identical planar 
3D uni-axial metamaterials that have only electric anisotropy, shows that the interaction is always attractive. 

A key issue here is the realization that it is incorrect to use the above Drude-Lorentz expressions when computing 


dispersion interactions (Rosa et al.. 2008bI. Indeed, although Eqs. (18) are valid close to a metamaterial resonance. 
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FIG. 10 (Color Online) Top: Normalized Casimir force as a function of separation between a gold semi-space and a 2D 

metasurface made of close-packed (lattice constant a = 11.24/im) LiTaOa spheres of radius r = 5.62/rm, calculated via scattering 

^2 _^2 

theory. The optical response of LiTaOs is described by a single-resonance Drude-Lorentz model, e(tu) = £oo(l + ^2 

where too = 13 . 4 , the transverse and longitudinal optical phonon frequencies are wt = 26.7 x 10 ^^ rad/sec and lul = 46.9 x 10 ^^ 
rad/sec, and the dissipation parameter is 7 = 0.94 x 10 ^^ rad/sec. The gold Drude parameters are taken as hQp = 3.71 eV and 
= 20 . Bottom: effective permittivity and permeability of the close-packed LiTaOs spheres as a function of (a) real and 


(b) imaginary frequencies as calculated by the Maxwell-Garnett effective medium theory. Figures taken from (Yannopapas and 


Vitanov 


20091. 


they do not hold in a broad-band frequency range. In particular, calculations based on Maxwell’s equations in a long 


wavelength approximation for SRRs result in a slightly different form for the effective magnetic permeability (Pendry 


et al. 1999) 


Msrr(w) = 1 - 




-b 


(19) 


where the filling factor / < 1 is a geometry dependent parameter. The crucial di fference betw een Eqs. (18) and 
(19) is the factor in the numerator of the latter, a consequence of Faraday’s law (Rosa 2009). Although close to 
the resonance both expressions give almost identical behaviors, they differ in the low-frequency limit: /reff(*w) > 1 
while < 1- The fact that all passive materials have e(i(jj) > 1 implies that Casimir repulsion is impossible 

for any magnetic metamaterial made of metals and dielectrics (Rosa et al.. 2008a). This conclusion is confirmed 


by scattering theory calculations, that do not rely on any effective medium or homogenization approximations. For 


example, in (Yannopapas and Vitanov 2009) the Casimir force was computed exactly for 2D metasurfaces made of a 


square close-packed array of non-magnetic microspheres of LiTaOs (an ionic material) or of CuCl (a semiconductor). 
Although the systems are magnetically active in the infrared and optical regimes, the force between finite slabs of 
these materials and metallic slabs is attractive since the effective electric permittivity at imaginary frequencies is 


larger than the magnetic permeability. In Fig. 10 we show the Casimir force (normalized with respect to the ideal 
zero-temperature Casimir force Fq) between a gold plate and a 2D LiTaOs metasurface together with the effective 
permittivity and permeabilities of a close-packed LiTaOs crystal. The results confirm that the Casimir interaction 
is attractive in magnetic metamaterials made of non-magnetic meta-atoms. In contrast, intrinsically magnetic meta¬ 
atoms could potentially lead to Casimir repulsion. Naturally occurring ferromagnets do not show magnetic response 
in the infrared and optical regimes, as needed by the Casimir effect, but small magnetic nanoparticles (e.g. few 
nanometer-sized Ni spheres) become super-paramagnetic in the infrared. A realization of the original idea for Casimir 
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repulsion by Boyer was then proposed based on a metasurface made of such intrinsically magnetic nano particles 


(Yannopapas and Vitanov, 2009). 


Chiral metamaterials made of metallic and dielectric meta-atoms were also proposed as candidates for Casimir 


repulsion (Zhao et al.. 20091. When described by an effective medium theory, such systems posses an effective magneto¬ 
electric response that modifies the standard constitutive relations in Maxwell’s equations as D = cgeE-|- zk^H/c and 
B = Close to a resonance, the magneto-dielectric coefficient can be modeled as Km{uj) = 

For such materials the reflection matrix is no longer diagonal and there is polarization mixing. 




Repulsive Casimir forces and stable nanolevitation was predicted for strong chirality (large values of uj^^lojK^r) (Zhao 
et al ., 20091. However, these results were shown to be incompatible with the passivity and causal response of the 
materials (Silveirinha and Maslovski 2010), which implies that the condition Im[e(a;)]Im[^(tLi)] — (Im[Km(a;)])^ > 0 
must be satisfied. This relation imposes a limit of the strength of the imaginary part of Km-, and results in an 
attractive Casimir force between chiral metamaterials made of metallic/dielectric meta-atoms for any physical values 
of the magneto-electric coupling ([Silveirinha 


2010). These theoretical arguments were also confirmed by full-wave 


simulations of chiral metamaterial structures (McCauley et al.\ 2010b|, and it was shown that microstructure effects 


(i.e. proximity forces and anisotropy) dominate the Casimir force for separations where chirality was predicted to 


have a strong influence. Still, chiral metamaterials may offer a way to strongly reduce the Casimir force (Zhao et al. 


2010 ). 


B. Photonic crystals 


Photonic crystal are man-made electromagnetic structures that, unlike metamaterials, have unit cell sizes on the 
order of the wavelength of light. The most important property of photonic crystals made of low-loss dielectric 


periodic structures occurs when the wavelength is about twice their period (Joannopoulos et al.. 2008). Many exotic 


phenomena are found, including the appearance of photonic band gaps preventing light from propagating in certain 
directions with specified frequencies, the localization of electromagnetic modes at defects, and the existence of surface 
states that bound light to the surface for modes below the light line. Photonic crystals were co-discovered in 1987 
by Eli Yablonovitch, who proved that spontaneous emission is forbidden when a three-dimensional periodic structure 
has an electromagnetic bandgap which overlaps with an electronic band edge (Yablonovitch 1987), and by Sajeev 


John, who showed that strong localization of photons could take place in disordered dielectric superlattices (John 


1987). The simplest possible photonic crystal, a ID multilayered stack made of materials of alternating dielectric 


constants, had been already investigated more than a century ago by Lord Rayleigh. Today photonic crystals come 


1991)), periodic 


in different fashions, including complex 3D structures (e.g., the Yablonovite (Yablonovitch et al. 
dielectric waveguides, and photonic-crystal slabs and fibers. 

Photonic crystals offer great flexibility in designing atomic traps close to surfaces at sub-micrometer distances 
allowing the integration of nanophotonics and atomic physics with a host of exciting quantum technologies. Trapping 
atoms near surfaces is determined by the Casimir-Polder force. However, in analogy to Earnshaw’s theorem of 
electrostatics, there are no stable Casimir-Polder (or Casimir) equilibria positions for any arrangements of non¬ 
magnetic systems, provided the electric permittivities of all objects are higher or lower than that of the medium 
in-between them (Rahi et al. 2010| ). For example, there is no stable equilibrium position for a ground-state atom 
above a metallic/dielectric structure]^ Fortunately, no such constraints exist for excited state atoms, or when the 
trapping potential energy is the superposition of the Casimir-Polder interaction and an external optical trapping field. 
The Casimir-Polder interaction can be calculated for an atom in state I with polarizability ai(uj) considering Eq. 
for a rarefield dielectric, as shown in (Dzyaloshinskii et al. 1961 Lifshitz and Pitaevskii 1980 Milonni[ 1993| ). 
Typically, the arising Green’s function is solved via computational FDTD techniques (to be reviewed in next section). 
It has been shown that the Casimir-Polder force between a ground-state atom and a ID dielectric grating can trap 
atoms along the lateral directions of the dielectric surface (Contreras-Reyes et al.. 2010). However, there is no 


trapping along the directions parallel to the grating’s grooves. Fully stable traps in 3D can be obtained utilizing 


photonic crystals, such as ID periodic dielectric waveguides (Hung et al. 2013). These proposed structures support a 


guided mode suitable for atom trapping within a unit cell, as well as a second probe mode with strong atom-photon 
interactions. The combination of the light-shifts from a laser beam together with the Casimir-Polder force from the 


^ A corollary of this theorem is that there is no Casimir repulsion for any metallic/dielectric-based metamaterial treated in the effective- 
medium approximation. Hence, when applied to dispersive interactions, effective medium is a good approach only at separations larger 
than the unit cell dimensions of the metamaterial. At short distances, displacements of structured Casimir plates might lead to repulsion 
that, however, must be compatible with the absence of stable equilibria. 
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d(nm) 


FIG. 11 (Color Online) Top: A) Schematic of the experimental configuration used to measure the Casimir force between a 
gold-coated sphere and a gold nanostructured grating. One of the many nanostructures used in the experiment is shown (SEM 
images in B-D). The radius of the Au sphere is R — 151.7 fim. Bottom: Plane-grating pressure as a function of sphere-grating 
separation for a lamellar grating with period p — 250 nm, tooth-width w = 90 nm, and height /i = 216 nm. Experimental 
measurements of sphere-grating force gradient divided by 2'kR (dots with error bars), and plane-grating pressure computed with 
the proximity force approximation (dashed lines) and exactly using scattering theory (solid line). Figure taken from ([Intravaial 


et al. 20131. 


dielectric nanostructure results in a fully stable, 3D atomic trap. Aligning the photonic band gap edges with selected 
atomic transitions substantially enhances the atom-photon interactions, since the electromagnetic density of state 
diverges due to a van Hove singularity. These ideas have been recently implemented experimentally with a Cs atom 
trapped within a ID photonic crystal waveguide consisting of two parallel SiN nanobeams with sinusoidal corrugation 
(Goban et al. 2014). The measured rate of emission into the guided mode along the ID waveguide was Fid = 0.32r', 
where F' is the decay rate into all other channels. Such a high coupling rate is unprecedented in all current atom- 
photon interfaces, and paves the way for studying novel quantum transport and many-body phenomena in optics. 
Other works involving atom-surface dispersive interactions in close proximity to photonic crystals include resonant 
dipole-dipole energy transfer (Bay et al. 1997) and enhanced resonant forces (Incardone et al. 2014) between atoms 
with transition frequencies near the edge of the photonic bandgap, and strong localization of matter waves mediated 


by quantum vacuum fluctuations in disordered dielectric media (Moreno et al. 2010b). 


C. Plasmonic nanostructures 


Metallic nanostructures can support collective electromagnetic modes, such as surface plasmons (also known as 
surface plasmon polaritons), which can propagate along the surface, decay exponentially away from it, and have a 
characteristic frequency of the order of the plasma frequency (Maier 2007). Surface plasmons affect the Casimir inter¬ 


action in a non-trivial manner (Intravaia and Lambrecht 2005), and this point was also discussed for Dirac materials 


in Sec. III. When written in terms of real frequencies, the Lifshitz formula, Eq. for planar systems has a term 
arising from the propagative modes, which gives an attractive force at all distances. There is a second term associated 
with the evanescent hybrid plasmonic modes, which results in an attractive force at short distances (shorter than the 
plasma wavelength) and a repulsive one at longer distances. There is a subtle cancellation between the attractive 
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and repulsive terms at large separations, resulting in an always attractive force between planar metallic surfaces 
for all separations. This observation suggests that metallic nanostructures at scales below the plasma wavelength 
can potentially enhance the repulsive contribution due to plasmons and lead to a suppression of the Casimir force. 
Nanostructured metallic surfaces with tailored plasmonic dispersions have already impacted classical nanophotonics, 


with applications ranging from extraordinary light transmission (Ebbesen et al. 19981 to surface-enhanced Raman 


scattering (Nie and Emory 19971. Metallic structures with strong deviations from the planar geometry and possessing 


geometrical features on very small scales are also likely to give significant new insights into potential Casimir devices. 

In addition to plasmons associated with the metallic nature of the plates, there is another type of plasmonic 
excitations, the so-called spoof plasmons, that arise from geometry and exist even for perfectly reflecting surfaces. 


Pendry and co-workers (Garcia-Vidal et al. 2005 Pendry et al. 2004) proposed engineered dispersion by periodically 


nanostructuring surfaces by perforating perfect electrical conductors. The resulting surfaces support surface modes 
that have dispersion similar to real surface plasmons in metals, but with the effective plasma frequency determined by 
the geometric parameters of the perforation. Spoof plasmons are also present in nanostructures made of real metals, 


enhance the modal density of states, and modify the Casimir interaction in nanostructured metallic cavities (Davids 


et al. 2014). 


Besides computations of the interaction in complex systems, including nanostructured surfaces (Biischer and Emig 


2004 Davids et al. 2010 Guerout et al. 2013| Intravaia et qL[ |201^ [Lambrecht and Marachevsk^q 2008 Noto et al. 


20141, advances in Casimir force measurements have also been reported. However, the experimental progress has 


been limited due to difficulties associated with the reliable fabrication and the measurement of the force. Using an 
in situ imprinting technique, whereby the corrugation of a diffraction grating was imprinted onto a metallic sphere 
by mechanical pressure, the lateral Casimir force between two axis-aligned corrugated surfaces was measured as a 
function of their phase shift (Chen et al. 2002), and the normal Casimir force between them was also measured as a 


function of the angle between their corrugation axes (jBanishev et al. 2013c). Nanostructured lamellar gratings made 


of highly-doped Si have been used to measure the Casimir interaction with a metallic sphere (Chan et al. 2008), 


with conclusive evidence of the strong geometry dependency and non-additivity of the force. More recently, a strong 


Casimir force reduction through metallic surface nanostructuring has been reported (Intravaia et al. 2013). In Fig. 


11 the experimental setup is shown, consisting of a plasmonic nanostructure in front of a metallic sphere attached 
to a MEMS oscillator. A deep metallic lamellar grating with sub-100 nm features strongly suppressed the Casimir 
force, and for large inter-surface separations reduced it beyond what would be expected by any existing theoretical 
prediction. Existing state-of-the-art theoretical modeling, based on the proximity force approximation for treating 
the curvature of the large-radius sphere (i? = 151.7/rm, much larger than any geometrical length scale in the system), 
and an exact ab initio scattering analysis of the resulting effective plane-grating geometry, did not reproduce the 
experimental findings. The development of a full numerical analysis of the sphere-grating problem, capable of dealing 


2013) with plasmonic nanostructures, 


with the disparate length scales present in the experiment (Intravaia et al. 
remains an open problem. 

Nanostructured surfaces have also been used in studies of atom-surface dispersion interactions. Casimir-Polder forces 
between a single atom or a Bose-Einstein condensate above a grating have been measured using different methods 


(Grisenti et al. 1999 Oberst et al. 2005 Pasquini et al. 2006 Perreault and Cronin, 2005 Zhao et al. 2008), and 


the near- and field-field scaling laws of the Casimir-Polder potential were verified. Theoretical proposals have also 
been put forward to measure the Casimir-Polder potential at corrugated surfaces with Bose-Einstein condensates 
(Messina et al. 2009 Moreno et al. 2010a). In a series of recent experiments, ultracold atoms were utilized to 


survey the potential landscape of plasmonically tailored nanostructures. In (Stehle et al. 2011) a Rb Bose-Einstein 


condensate was accelerated towards Au plasmonic microstructures whose plasmons were excited by external laser 


fields in a Kretschmann configuration (Fig. 12). A blue-detuned laser beam generates an evanescent optical field that 
repels the atoms from the surface, while the atom-grating Casimir-Polder interaction produces an attractive potential. 
This combination results in a potential barrier that can be mapped by classical or quantum reflection measurements. 
Diffraction measurements of Bose-Einstein condensates from metallic nanogratings have allowed to locally probe the 
Casimir-Polder potential ( ]Bender et al. 2014), revealing information about its landscape (Fig. [l^ in agreement with 
theoretical calculations based on the scattering approach to atom-grating Casimir interactions. 


V. NON TRIVIAL BOUNDARY CONDITIONS 


Casimir interactions are fundamentally changed in the presence of macroscopic objects due to the shapes of bound¬ 


aries and interfaces, which lead to complex and highly non-additive wave effects (Buhmann 2012b Dalvit et al ., 2011a 


Reid et al. 2013a Rodriguez et al.. 2011a 2014). Understanding the ways in which non-trivial shapes and boundary 
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FIG. 12 (Color Online) Top: Schematic configuration to measure the Casimir-Polder potential probed by a Bose-Einstein 
condensate diffracted from a plasmonic nanostructured being excited by an external laser field in the Kretschmann configuration. 
Bottom: Atom-grating potential landscape arising from the combination of a repulsive evanescent-wave potential and the 
Casimir-Polder attraction. The external laser field has a power of P = 211 mW. The evanescent field from the grating 
modulates the repulsion, and the attractive Casimir-Polder potential is the strongest on top of the gold stripes. At a distance 
of 200 nm, the potential is laterally modulated with an amplitude of AE/ks = 14/iK. Figures taken from ([Bender et ol.[[2014| 


Stehle et al. 20111. 


conditions affect the force has not only shed light on various ways to design forces used to combat unwanted Casimir 
effects in NEMS and MEMS, but continues to reveal regimes and situations where the often-employed proximity force 
approximation (PEA) and pairwise summation (PWS) approximation fail dramatically. Such structures can also lead 
to forces that differ significantly from the attractive, monotonically decaying force laws associated with planar bodies 
and/or dilute, atomic media. 

Early studies of Casimir forces focused on simple geometries, e.g. planar bodies and generalizations thereof, by 
employing sum-over-mode formulations where the zero-point energy of electromagnetic fields (field fluctuations) rather 
than dipolar interactions (charge fluctuations) were summed (Casimir 1948| Milonni 1993). The equivalence of these 


two perspectives comes from the fluctuation-dissipation theorem, relating the properties (amplitude and correlations) 


1984 

Lifshitz 

1956 

Lifshitz and Pitaevskii 

1980 


1980). Ultimately, the connection between current and field fluctuations 


arises from the well-known dyadic electromagnetic Green’s function (Jackson 1998): 


Gij{r, r/w) = I [V X V X —e(r,a;)cLi^] ^ ej(5(r — r')| 


( 20 ) 


where is the unit vector. 


The connection to sum-over-mode formulas arises from the trace of the Green’s func¬ 
tion being related to the electromagnetic density of states p{ijj) = /^ Tr Im Gj, (r, r, cc), which when integrated 


p{ijj) = J dujp{u}) leads to the famous £ = ^ formula (Gerlach 


1971 


van Kampen et al. 


1968 


Rodriguez 
y metallic 

conductors, where Hermiticity leads to well-defined modes, it has also been extended to handle other situations of 


et al.\ 2007b I. Although this formulation was originally developed in special geometries involving perfect 


interest such as open structures and lossy dielectrics ( 

Davids et al. 

2010 

Genet et al. 

r 

>003 

Graham et al. 

2009 

Intravaia and Behunin 

2012 

Milton et al. 

2010b Mochan and Villarreal 

2006 

Van Enk 

1991 

)). Despite these gen- 


eralizations, the sum-over-mode approach poses practical challenges for computations in general structures due to the 


cumbersome task of having to compute all of the modes of the system (van Enk 1995 Ford 1993 Rodriguez et al. 


2007b I. 


Instead, more powerful applications of the fluctuation-dissipation theorem exist in which Green’s functions are 
directly employed to compute energy densities and stress tensors (momentum transport) rather than modal contri¬ 
butions to the energy, reducing the problem to a series of classical scattering calculations: scattered fields due to 
known incident fields/sources. This latter viewpoint was originally employed by Lifshitz and others to calculate forces 
between planar dielectrics bodies ( [Dzyaloshinskh et al. 1961 Lifshitz and Pitaevskii [1980 ), and it turns out to be 
much more useful when dealing with complex geometries. The advantage comes from the fact that the Green’s func¬ 
tion does not need to be obtained analytically, as was done for planar bodies, but it can be routinely and efficiently 
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FIG. 13 (Color Online) Schematic illustration of numerical methods recently employed to compute Casimir interactions in 
complex geometries: (a) scattering methods where the field nnknowns are either incident/outgoing propagating waves (left) 
or vector currents J(x) defined on the surfaces of the bodies (right), and the resulting energies are given by Eq. 21 (b) 


stress-tensor methods in which the force is obtained by integrating the thermodynamic Maxwell stress tensor over a surface 
surrounding one of the bodies utilizing the fluctuation-dissipation theorem. 


computed numerically via classical electromagnetism. These ideas lie at the center of recently developed general- 
purpose techniques, schematically shown in Fig. |13[ for computing forces in complex structures that boil down to 


a series of classical scattering calculations of Green’s functions (McCauley et al. 

2010a Pasquali and Maggs 

2008 

2009 Rodriguez et al. 

2007bl|2009 Xiong and Chew 2009 Xiong et al. 2010) or related quantities such as scattering 

matrices (. 

Atkins et al. 

2013 iBiischer and Emig 2004 Emig 

2003 

Emig et al. 

2007 Gies and Klingmuller 

2006b 

Gies et al. 

2003 Kenneth and Klich 2008 Lambrecht et al. 

2006 

Milton et al. 

2010al Rahi et al. 2009 Reic 

et al. 

2011 2009 2013b). 


Despite their relative infancy, these methods have already led to a number of interesting predictions of unusual 


Gasimir forces in a wide range of structures, including spheres (E 

Imig 

2008 Gies and Klingmuller 

|2006a 

Maia Neto 

et al. 2008 

, cylinders 

(Emig et al. 2006), cones ( 

Maghrebi et al. 

2011 

), waveguides (Pernice et al 

( 2010 

Rahi et al. 

2008 

Rodriguez et al. 

2007a Rodriguez-Lopez et al. 2009 Zaheer et al. 2007), and patterned surfaces (1 

Biischer and 

Emigl 

2004 

Davids et al.. 2010 Emig et al. 2003 

Guerout et al. 

2013 

Lambrecht and Marachevsky, 2008 Rodrigues 

et aL\ 

2006 

Rodriguez et aL\ 2008a), among others (Broer et al. 

2012 

Rodriguez et al.\ 2010a 2008b, 2010b). Here 


we provide a concise but inclusive exposition of the main techniques employed in state-of-the-art calculations along 
with discussions of their suitability to different kinds of problems, all the while focusing on representative results that 
reveal the highly non-additive character of Casimir forces. 


A. Scattering methods 


A sophisticated and powerful set of techniques for calculating Casimir forces are scattering methods (Fig. (Ha)). 
While these approaches come in a variety of flavors, they often rely on formulations that exploit connections between 


the electromagnetic density of stat 
path-integral representations of the 

es (Emig 2003 Emig et al. 2006 

Gies et al. 2003 Kenneth and Klich 

|2008|) or 

electromagnetic energy (Emig et al 

2007 Kardar and Golestanian 

1999 

Kenneth 

and Klich, 2008 Lambrecht et al. 

2006 

Milton et al. 2010a Rahi et al. 2009 Reid et al. 2011 2009 

), anc 

[ classical 

scattering matrices (Dalvit et al. 

2011a 

). Regardless of the chosen starting point, the Casimir energy turns out to 


be given in the form: 


ft 

£ = -/ dw logdet(MMoi), (21) 

27r Jo 

where the matrix M is 

M(11) ... 

M(22) . . . 



whose diagonal blocks 


are precisely the scattering matrices of isolated bodies and whose off-diagonal blocks 


(“/3) encapsulate interactions and scattering among the bodies (Lambrecht et al. 


2006 


Rahi et al. 


2009). Mul- 
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ensures that the divergent, self-interaction energy of the bodies in isolation 


tiplication by the inverse matrix 
(separations d —> oo) is subtracted, leaving behind a finite quantity (Rahi et al. 


n[2c 
• & 


2009). 


Although Eq. 21 may appear largely unrelated to the Lifshitz formula (Eq. It^, the connection between the two 
becomes apparent when considering extended structures. Specifically, given two semi-infinite, periodic bodies the 
formula can be written in the more familiar form, £ = dio logDet (l — , where are the 

reflection matrices of each individual half-space and are translation matrices that describe wave propagation 


between them (Lambrecht et al. 2006 Lambrecht and Marachevsky 2008 Rahi et al.. 2009). For planar bodies, as 


discussed in Sec. III.A, the scattering matrices can be expressed in a Fourier basis and the above expression reduces 
to the Lifshitz formula (Lifshitz 1956), originally obtained via direct evaluation of the Maxwell stress tensor. 


It is also possible to derive a slightly different scattering formula, known as the TGTG formula (Kenneth and Klich 


2008 Klich and Kenneth, 2009), in which the energy between two arbitrary bodies is expressed as 


£ = 


dn logdet (l - 


( 22 ) 


where are the T operators appearing in the Lippmann-Schwinger equation (related to the scattering matrices of 
individual bodies), and Gq“^^ are the homogeneous Green’s functions of the intervening medium, describing the wave 
propagation. Note that even though these formulations may appear to be completely divorced from the original picture 
of dipole fluctuations, the fact that the energy is described by the scattering properties of the bodies is not surprising. 
In particular, as discussed further below, at equilibrium it is possible to describe the statistics of Held fluctuations 


independently of the corresponding current sources of the fluctuations (Eckhardt 1984 Lifshitz and Pitaevskii, 1980). 


Intuitively, one can consider Casimir interactions as arising from the scattering and momentum-exchange of vacuum 
electromagnetic fields originating from radiating sources infinitely far away (rather than within the bodies) and that 
ultimately end up equilibrating as they get scattered, absorbed, and re-emitted by the bodies. 

Equationj^was originally exploited to study forces between highly symmetric structures, e.g. spheres and cylinders, 
where the corresponding propagators, scattering, and translation matrices can be expanded in terms of convenient. 


de-localized free-wave solutions of the Helmholtz equation (Kenneth and Klich 2008 

Lambrecht and Marachevsky 

2008 Maghrebi et al. 2011 Milton et al. 

2010a [Rahi et al. 

2009), as illustrated in Fig 

13[a). The resulting spectral 

methods (Balian and Duplantier 1978 ' 

Jalvit et al. 2006 

Emig et al. 

2007 Kenneth and Klich, 2008 

Lambrecht 

et al. 2006 Lambrecht and Marachevsky 2008 Mazzitelli et al. 2006 

Milton et al. 

2008| [Rahi et al. 

2009) are 


advantageous in a number of ways: First, they yield analytical results that offer insight into the properties of the 


Gasimir force at asymptotically large separations (Rahi et al. 2009) or under assumptions of dilute media (Bitbol 


et al. 2013 Golestanian 2009 Milton et aZ. 2008 Milton and Wagner 2008). Second, the trace operations for smooth 


and high-symmetry structures can be efficiently implemented due to the very high-order and possibly even exponential 


convergence of the basis expansions (Boyd 2001 Dalvit et al. 2011a). Finally, since the energy expressions involve 


simple products of scattering matrices having well-studied properties, this formulation is well-suited for establishing 
general constraints on the signs and magnitudes of forces under various circumstances. Of particular importance is the 


recent demonstration that the force between any two mirror-symmetric bodies must always be attractive (Kenneth 


and Klich 2006 Kenneth et al. 2002), resolving a long-standing question about the sign of the internal pressure or 


self-force on a perfectly metallic, isolated sphere (the limit of two opposing hemispheres) (]Bordag et al. 2001| Boyer 


1968[ Brevick and Einevoll 1987 Milton et al.. 1978). Similarly, recent works have shown that stable suspensions 


(local equilibria) between vacuum-separated, non-magnetic bodies are generally impossible (Lambrecht et al. 1997 

2oTo| . 


Rahi et al. 


For more complicated bodies lacking special symmetries, involving sharp corners, or where non-uniform spatial 
resolution is desired, it is advantageous to employ localized basis functions. More commonly, the unknowns are 
defined on a generic mesh or grid and the resulting equations are solved numerically, examples of which are the finite 


difference (Taflove and Hagness 2000), finite element (Jin 2002), and boundary-element (Bonnet, 1999 Chew et al. 


2001) methods. The latter category are closely related to scattering methods (Reid et al.. 2011 2009 Rodriguez et al. 


2007b Xiong et al. 2010); in the surface-integral equation formulation of electromagnetic scattering, the scattering 
unknowns are fictitious electric and magnetic currents defined on the surfaces of the bodies, illustrated in Fig, [l^a), 
and expanded in terms of an arbitrary basis of surface vector fields (Chew et al. 2001). The connection to scattering 


problems comes from the fact that incident and scattered fields are related to the current unknowns via homogeneous 
Green’s functions (analytically known); not surprisingly, this formulation leads to a similar trace expression for the 
Casimir energy given in Eq. except that the elements of M consist of overlap integrals among the various surface 
basis functions. A powerful implementation of this approach is the boundary-element method (BEM), where the 
current unknowns are expanded in terms of localized basis functions (typically, low-degree polynomials) defined on 
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the elements of some discretized surface, as illustrated in Fig. 13’a). As a result the 


matrices turn out to be none 
Such a 


20011 . 


other than the well-studied BEM matrices that arise in classical scattering calculations (Chew et al. 
formulation allows straightforward adaptations of sophisticated BEM codes, including recently published, free and 
widely available software packages (Reid, 2012). Like most numerical methods, the BEM method can handle a wide 
range of structures, including interleaved bodies with corners, and enables non-uniform resolutions to be employed as 
needed. 


B. Stress tensor methods 


Although originally conceived as a semi-analytical method for computing forces in planar bodies (Jaekel and 
Reynaud 1991[ Klimchitskaya et al. 2000 Tomas, 2002 Zhou and Spruch 1995), leading to the famous Lif- 


shitz formula (Sec. Ill.A), the stress tensor approach can also be straightforwardly adapted for numerical com¬ 
putations (McCauley et al. 2010a Rodriguez et al. 2007b 2008a 2009[ 2008b[ Xiong and Chew 2009) since 
it relies on repeated calculations of Green’s functions. In this formulation (schematically shown in Fig. [Id^b)), 
the Casimir force on an object is expressed as an integral of the thermodynamic, Maxwell stress tensor (Ty ) = 
e (^{EiEj) — i -I- (^{HiHj) — ^ over an arbitrary surface S surrounding the object. 


F = / dK / / (T) • dS. 

Jo J Js 


(23) 


Similar to scattering methods, here the picture of fluctuating dipoles is masked by an equivalent scattering problem 
involving fields rather than fluctuating volume currents, whereby the correlation functions (EiEj), {HiElj) ~ Gj^, a 


consequence of the fact that at equilibrium currents and field fluctuations become thermodynamically equivalent (Eck- 
hardt |1984 ). 


Beyond special-symmetry structures where the Green’s functions can be expanded in a convenient spectral ba¬ 


sis (Jaekel and Reynaud 1991 Klimchitskaya et al. 2000 Tomas 2002 Zhou and Spruch 1995), recent implemen¬ 


tations of the stress-tensor method for arbitrary geometries exploit general-purpose techniques, such as the finite- 
difference method illustrated in Fig. mb), where space is divided into a uniform grid of finite resolution, and the 


resulting matrix equations for the Green’s functions are solved numerically (Anderson et al. 1999 Strikwerda, 1989 
Taflove and Hagness, 2000). Early implementations include both finite-difference frequency-domain (jRodriguez et oL 


2007b Xiong and Ghew 2009) and time-domain (McCauley et al. 2010a Rodriguez et al.. 2009) methods. 


Because Casimir forces involve broad bandwidth fluctuations time domain methods are advantageous in that 


Gij(r,r',w) at all frequencies can be computed at once via Fourier transforms (Taflove and Hagness, 2000). While 


the finite difference stress-tensor method does not offer the efficiency and sophistication of other formulations and 


discretization schemes, such as the BEM fluctuating-surface current method (Reid et al. 2011), they are compensated 


by their flexibility and generality. For instance, they are extremely simple to implement (leading to many free and 


easy-to-use numerical packages (Oskooi et al.. 2010), can handle many different kinds of boundary conditions and 


materials (including anisotropic and even nonlinear dielectrics), and are well understood. A BEM implementation 


of the stress-tensor method was also first suggested in Ref. (Rodriguez et al. 2007b) and subsequently implemented 


by Ref. (Xiong and Chew 2009), although for small problems the trace formulas above provide a simpler and more 


efficient alternative since they do not require repeated integration over surfaces and involve only products of BEM 
matrices. On the other hand, the stress tensor method offers computational advantages for large problems since 
it involves repeated evaluation of Green’s functions, or matrix-vector products, making it an ideal candidate for 


applications of fast-solver (iterative) techniques (Chew et al. 1997). 


C. Casimir interactions in complex geometries 


While PFA and PWS approximations provide simple, quickly solvable, and intuitive expressions for forces in ar¬ 
bitrary geometries, they are uncontrolled when pushed beyond their limits of validity and have been shown to fail 


(even qualitatively) in the simplest of structures (Bitbol et al. 2013 Bordag 2006 Dalvit et al.. 2011a Gies and 


Klingmuller 2006a[ Rodriguez et al. 2011a). Increased demand for experimental guidance has stimulated recent 


efforts in quantifying the validity and accuracy of PFA. 

Although PFA is technically only applicable in geometries with smooth, large-curvature objects and small separa¬ 


tions, it has nevertheless been heuristically applied in the past to study a wide range of other situations (Lambrecht 
et al.. 2008 Rodriguez et al. 2014). For instance, a number of recent works have employed scattering methods 
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FIG. 14 (Color Online) Results illustrating the range of validity of PFA and PWS approximations in the sphere-plate geometry, 
involving a sphere of radius R separated from a semi-infinite plate by a surface-surface distance d. The exact Casimir energy is 
computed by application of the scattering method, (a) Ratio of the exact and PFA energies as a function of d/R for perfectly 
metallic conductors ( [Bimonte et Numerical data (dots) are compared to the first correction of the PFA (dashed 

lines), described in the text, and also to a fit performed via Fade approximants (solid curves), in which the energy ratio is 
given by + + TSo S' (^) of f ^6 PWS and exact Casimir energies as a function of the static 

permittivities of the sphere and plate, for various ratios d/R (blue curves) (Bitbol et al. 20131. The red and black curves 


represent equivalent results for the plate-plate and sphere-plate geometries in the limit of infinite separations. 


to investigate extensions of PFA in the sphere-plate geometry at large separations in the idealized limit of perfect 
conductors (Bimonte et al. 2012 Fosco et al.. 20121, where the PFA energy takes on the closed-form expression 
SpFA = 1440 ^ ; where R denotes the sphere radius. The plotted ratio of energies in Fig. 14 'a) shows that PFA increas¬ 
ingly overestimates the energy as d —> oo, which is to be expected si nce at d/R^ 1, the interacti on approaches that 

Higher-order per- 


of a dipole above a plate exhibiting a significantly faster decay ~ Jg (Buhmann and Welsch 


2007 


turbative PFA corrections for large curvatures d/R 1 have also been obtained ([Simonte et a/.[|2012 Bordag, 2006 


Fosco et al. 2012 Lambrecht et al. 20081 (Fig. [l4[a)). Techniques based on Fade approximants, which constrain 


the force at short and large separations using gradient and multipole expansions (Fig. [I4|(a)), have also culminated 


in analytical expressions (Bimonte et al.. 2012), leading to the hope that similar methods can be applied to more 
complex geometries. The situation is more complicated in cases involving realistic metals and finite temperatures, as 
illustrated by recent predictions of geothermal effects involving nontrivial interplay between geometry, materials, and 


temperature in the sphere-plate geometry (Canaguier-Durand et al. 
2010a|b| ). 


2010 Maia Neto et al.. 2008 Weber and Gies 


The PWS approximation, applicable in the limit of large separations and dilute media, relies on dividing the object 


into small elements (” atoms”) and summing the corresponding vdW and Casimir-Polder interactions (Bergstrom 1997 


Golestanian 

2009| Milton et al. 

to 

o 

o 

00 

Veble and Podgornik 

2007). The presence of multiple scattering, otherwise 

absent in the limit of weak coupling or dilute media ( 

Milton et al. 

2008 

), has long been known to signihcantly modify 

the underlying two-body force laws (Axilrod and Teller 1943 

). Despite these shortcomings, PWS approximations have 


been recently applied to numerically approximate interactions in complex geometries (Sedmik et al ., 2007 Tajmar 


20041, especially in the field of microfluidics (Parsegian 2006 Stone and Kim 2001). While PWS approximations are 


strictly applicable in the limit of dilute media, recently they were shown to lead to larger errors in the experimentally 
relevant case of dielectric materials (Bitbol et al. 2013). This situation is illustrated in Fig. [T4|(b), which shows 
that PWS underestimates the energy in the perfect-metal limit (e —)■ —oo) by « 20%, is exact in the dilute limit of 
e —>■ 1, and is (surprisingly) most inaccurate at intermediate e ^ 10 where it overestimates the energy by roughly 
60%. Such counter-intuitive results shed light on the complexities associated with dilute approximations, since a 
heuristic argument based on the screening of fields in materials with large dielectric contrasts would predict strictly 
monotonically increasing deviations. By examining interactions between compact objects at asymptotically large 
separations, it is also possible to obtain perturbative corrections to Casimir-Polder forces ([Balian and Duplantier 


1977 


et al. 


EmigI |2008[ |Emig et~aL[ |2006[ |Golestanian[ |2000[ |2009[ [Milton and Wagiierl |2008| [Rahi et aL[|2009| Stedmanj 


2014). Formal derivations of PWS approximations in the limit of dilute media as well as perturbative corrections 


applicable in systems with larger index contrasts have also been developed (Golestanian, 2009 Milton et al. 2008 


Rodriguez-Lopez 2009). 


At intermediate separations that are on the order of the sizes of the objects and for realistic materials, neither PFA 
nor PWS, nor perturbative corrections thereof can accurately predict the behavior of the Casimir force. However, it is 
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FIG. 15 (Color Online) Selected results of Casimir interactions involving compact bodies suspended above planar objects 
or interacting with other compact bodies, calculated via scattering matrix techniques in combination with spectral methods, 
(a) Casimir energy of a perfectly conducting, vertically oriented cone of semi-opening angle 6o suspended above a perfectly 

2011|. The PFA approximation is shown as the dashed line, (b) Casimir 


conducting plate by a fixed distance d {Maghrebi et i 


energy of a perfectly conducting wedge and plate, with one face at a Hxed angle (po + ^o, as a function of the full opening angle 
= 26o (Maghrebi et al. 2011| . The PFA prediction, in dashed lines, is compared to the exact calculation, indicated by solid 
circles, (c) Ratio of the exact Casimir and PFA energies of the perfectly conducting cylinder-plate structure shown in the inset, 
decomposed into both TE and TM contributions (Emig et al. 20061. (d) Casimir force between perfectly conducting metallic 


cubes (red squares) or spheres (green circles), as computed by the BEM method of Ref. (Reid et al. 2013cI, divided by the 


corresponding PFA forces. Results for spheres are compared to computations performed using scattering matrix methods (blue 
circles). 


precisely this regime that is most easily tackled by numerical methods. Application of scattering methods to the study 
of compact bodies interacting with planar objects have led to a number of interesting predictions, a select number 
of which are illustrated in Fig. |15| In geometries involving perfect-conductor bodies with special symmetries such as 
spheres, cones, wedges, or cylinders, scattering-matrix methods have been employed to obtain both numerical and 


semi-analytical results (Dalvit et aZ. 2011a Emig et aZ. 2009 2006 Maia Neto et aZ. 2008 Mazzitelli et aZ., 2006 Rahi 


et al. 2009). Other, more complicated shapes such as waveguides, disks, cubes, tetrahedral particles, and capsules are 


less amenable to spectral methods, but have nevertheless been studied using brute-force techniques (Atkins et al. 2013 


Reid et al. 

2009 

2013c 

Rodriguez et al. 


9o and a substrate vanishes logarithmically £ ~ ~~ I losOn t 0, a type of divergence that is characteristic of lines 

and other scale-invariant objects (Maghrebi et al. I [2^ . In contrast, the PFA energy is predicted to vanish linearly 
as 9q —>■ 0. For a tilted wedge, the PFA energy remains constant until the back surface of the wedge becomes visible 


to the plate, while exact results indicate smoothly varying angle dependence despite the screening effects (Maghrebi 
et al. 2011). Earlier calculations of forces between cylinders, spheres and ellipsoidal bodies and plates have also 


demonstrated unexpectedly weak decay rates and other interesting non-additive modihcations (Dalvit et al ., 2006 


et al. 2013c I. 


Emig et aZ. 2009 ,2007 2006 Maia Neto et aZ. 2008 Mazzitelli et aZ. 2006). For more complicated structures, such as 


the pair of cubes shown in Fig. 15 it is more convenient to employ brute-force techniques like the BEM method (Reid 


Unusual Casimir interactions in multi-body geometries have also been recently studied (Dalvit et al. 2011a Ro¬ 


driguez et al ., 2011a). For instance, application of numerical methods (first employing stress tensors (Rodriguez et al. 
2007a) and subsequently scattering matrices (Rahi et 'al] [2008 )) in a structure composed of two metallic co-planar 


waveguides suspended above adjacent metal sidewalls (Fig. |16[a)) reveal that the attractive Casimir force per unit 
length between the waveguides varies non-monotonically as a function of their separation from the sidewalls h. Large 
deviations from PFA can be explained from the fact that PFA is unable to accurately capture the competing effects 
of TE and TM fields at small and large h (Hertzberg et aZ.l|2007l|Rahi et aZ.l|2008l|Zaheer et aZ.l|2007|). Extensions of 


this geometry to situations involving finite rods, such as the cylindrically symmetric geometry of Fig. |l^ a) where the 
sidewalls are joined to form a cylindrical tube and described by either perfect-electric or perfect-magnetic boundary 
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FIG. 16 (Color Online) Selected results illustrating unusual Casimir effects from non-additive interactions in complex structures, 
(a) Casimir force between either translationally invariant waveguides (solid lines) or cylindrically symmetric rods (dotted 
lines), normalized by the corresponding PFA force, as a function of their separation from adjacent sidewalls a distance h 
apart (McCauley et al. 2010a Rodriguez et al. 2007aI. Configurations of either perfect-electric (blue) or perfect-magnetic 


(red) conductors are considered, (b) Casimir pressure between two patterned structures involving periodic arrays of cylinders 
embedded in a semi-infinite substrate, as a function of their surface-surface separation d, computed by application of both 
scattering and FDTD methods (McCauley et al. |2010a 2011b Rodriguez et al. 20091. Exact calculations (solid lines), PFA 
(dotted lines), and effective-medium theory (dashed lines) are also shown. (F) Casimir force between a small, anisotropic 


particle or an array of particles and a plate with a hole ( Levin et al.\ 2010 McCauley et al. 
for vacuum-separated metallic bodies. 


2011al, demonstrating repulsion 


conditions (McCauley et al.. 2010aI, demonstrate the importance of dimensionality and boundary conditions on the 


behavior of the force. Along similar lines, structures involving periodic arrays of finite cylinders on slabs (Fig. |16[ b)) 
reveal strong variations in the force depending on whether the arrays are aligned or crossed, even leading to changes in 
its sign at close separations when the system is immersed in a fluid (McCauley et al. 2010a 2011b Rodriguez et al.\ 


20091. Exact calculations are also compared to predictions based on PFA, showing significant, qualitative deviations. 
An effective medium theory description of the problem, in which the slabs are treated as homogeneous, anisotropic 
dielectrics, gives surprisingly accurate predictions down to separations of the order of the period. 


Objects with nontrivial geometry can also be utilized to obtain repulsive Casimir interactions in vacuum (Rodriguez 


et al. 2014). A proof-of-principle of the feasibility of repulsion in vacuum was recently demonstrated using BEM and 


FDTD numerics in a structure involving a small, elongated particle above a plate with a hole (Levin et al. 2010), shown 


schematically in Fig.[l6|(c). Due to constraints on the size of the particles and hole as well as on the lengthscales needed 
to observe these effects, the force in that geometry turns out to be too small (atto-Newtons) for current experimental 
detection (Levin et al.[|2010|. However, extensions to multiple particles (attached to a substrate) have demonstrated 


thousand-fold force enhancements without the need to change hole radii or lengthscales (McCauley et al. 2011aI, 


as illustrated in Fig. [lejc). Interestingly, the enhancement can be understood as arising not only from the presence 
of additional bodies, but from increased repulsion due to the larger polarizability of the particles as they interact 
with fringing fields near the edge of the plate (Eberlein and Zietal 2011 Milton et al. 2011 2012). One can also 


show that the interaction between a polarizable particle and a perfect-metal wedge or half-plate is repulsive (Milton 
et al. 2012), provided that the wedge is sufficiently sharp and that the particle is sufficiently anisotropic. Similar 


results should extend to vdW interactions on molecules and atoms near structured surfaces, but the main challenge 
in these systems is the need to attain a large degree of particle anisotropy. Recent calculations show that Rydberg 
atoms cannot achieve a high enough anisotropy (Ellingsen et al. 2010). Regardless of their current experimental 


observability and practical considerations, these recent theoretical predictions demonstrate that geometry can prove 
to be a powerful resource for shaping Casimir forces. 


VI. SOFT AND BIOLOGICAL MATERIALS 

Besides the prominent role of fluctuation-induced interactions in inorganic materials systems, vdW forces have 
many other interesting manifestations. The adhesion of the gecko, with no help from glues, suction or interlocking, is 
perhaps the most popular example for vdW interactions in biological and bio-related matter. Researchers have shown 
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experimentally that vdW interactions between the gecko spatular toes and hydrophobic surfaces in air are responsible 


for the gecko clinging to substrates (Autumn and Peattie 2002 Autumn et al. 2002 Lee 2014). Other experiments 


suggest that while geckos indeed use no glue, they do leave ’’footprints” of residue identified as phospholipids with 


phosphocholine head groups (Hsu et al. 2012). On the other hand, the contact surface between the gecko’s toes and 


the substrate is saturated with methylene moieties of the phospholipids and contains no water. This is furthermore 


consistent with predominantly hydrophobic surface of gecko setae (Stark et al. 2013), landing additional support to 


the mostly vdW origin of the adhesion force. Currently there is much gecko-inspired interest in constructing materials 
with similar adhesion properties (Jeonga and Suh 2009). Single strand vertical arrays of cylindrical pillars produced 
by electron-beam lithography and etched into an array of vertical round shaped pillars are expected to show behavior 
similar to gecko toes. Dry adhesives for robotic applications based upon the characteristics of vertical and angled 
flaps from polydimethylsiloxane (PDMS) are also being considered for applications (Yu et al. 2011). 


A. The importance of aqueous solvent 


The most important defining characteristics of vdW interactions in soft- and bio-matter comes from the presence 
of a solvent, i.e., water ( Israelachvili| 1991| ). The interaction between two substrates with dielectric function s{u}) 
separated by a water layer with ew(uj) can be described by the standard Lifshitz formula (Eqj^ ( [Bordag et al.\ 
2009b ). The fact that typically £:(0) <C £^,(0) makes the n = 0 Matsubara term quite important, and can account for 
about 50 % or more of the total value of the Hamaker coefficient (Ninham and Parsegian[ |1970b[). For lipid-water 


systems (Kollmitzer et al. 2015) retardation effects are not important even at very large distances (Ninham and 


Parsegian, 1970a Parsegian 2006). The significant thermal effects from the n = 0 term show that results for vdW 


interactions in standard condensed media cannot be simply transcribed into the soft-matter context. Also, apart from 
its large static dielectric constant, the full dielectric spectrum of water leads to non-monotonic features in the vdW 
interaction between ice and water vapor (Elbaum and Schick 1991[) or hydrocarbon films (]Bar-Ziv and Safran 1993) 


across a liquid water layer in the retardation regime (Wilen et al. 1995). 


Solvent effects are also important when their dispersion properties are in a certain relation with those of interacting 


materials. Recent work by Capasso et al. (Munday et al. 2009 Munday and Capasso 2007), as well as previous work 
by various authors (Feiler et al. 2008 Lee and Sigmund, 2001, 2002 Meurk et al. 1997 Milling et a7[[l9% ), made 


it clear that for specific asymmetric interaction geometries, a solvent whose dielectric permittivity £m(w) is between 
those of the interacting bodies 1 and 2, £i(a;) > £m(w) > £ 2 (w), can create repulsive vdW interactions. Though in 
principle this solvent-mediated Casimir-Lifshitz levitation has been known since the appearance of the Lifshitz theory 
(Dzyaloshinskii et al. 1961), it has not been used to effectively control the sign of the vdW interaction. Solvent 


mixtures with low molecular weight solutes such as glucose and sucrose also affect the dielectric properties of the 


solution and can thus modify the vdW interactions (Neveau et al. 1977). Solvent-like effects could be important also 


for two graphene sheets separated by atomic hydrogen gas, with one sheet adsorbed on a Si 02 substrate, while the 
other is freestanding ( Bostrom and Sernelius[ 2012), as discussed in Sec. IIIC. 

Electrolyte screening is also a defining feature for bio-matter in aqueous environments. The presence of salt ions 
screens the n = 0 Matsubara frequency in the Hamaker coefficient. The existence of this screening is connected with 
the fact that the n = 0 Matsubara term actually corresponds to the classical partition function of the system, and for 
confined Coulomb fluids, such as inhomogeneous electrolytes, can lead to a thoroughly different form of the n = 0 term 
in the full Matusbara sum of the Lifshitz theory (Naji et al ., 2013 Podgornik and Zeks[ 1988), as discussed below. The 
screening of vdW interactions in electrolyte solutions can be derived in a variety of ways, most simply by replacing the 


Laplace equation with the linearized Debye-Huckel equation (Israelachvili 1991 Parsegian, 2006). In this approach. 


the free ions present in the aqueous solution are taken into account just like in the case of bad conductors (Pitaevskii 


2008 ), where the number of charge carriers is small and obeys the Boltzmann statistics. The zero frequency Matsubara 
term in Eq. [^then leads to an approximate n = 0 Hamaker coefficient 'Ho{d) = {3/4:)kBT{l + 2KQd)e~^'^°‘^ (Parsegian 


20061). This expression is obviously screened with twice the Debye screening length Kq ^ = SiriBno, where Ib — 0.7 


nm is the Bjerrum thermal length and no is the bulk salt concentration. 

In many bio-systems, the vdW free energy is typically cast into the form of the Hamaker-type approximation 
T{d,T) = with the separation-dependent Hamaker coefficient that can be calculated exactly via the Lifshitz 

formalism, when accurate experimental data for the dielectric properties are available. For example, the Hamaker 
coefficients for lignin and glucomannan interacting with cellulose, titania and calcium carbonate in vacuum, water 
and hexane, are found within a relatively narrow range of ~ 35 — 58 zJ for intervening vacuum and ~ 8 — 17 zJ for an 
intervening a queous medium ([Hollertz et al. 2013), with the dielectric response properties extracted via spectroscopic 
ellipsometry (Bergstrom et al. 1999). The Hamaker coefficients for the interactions of the wood components with 
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FIG. 17 (Color Online) Multilamellar array of lipid membranes, each composed of a hydrocarbon core with two surface 
hydrophillic headgroup layers. Left: a realistic presentation with fluctuating positions of the membranes (for details see 
(Kollmitzer et al\ 20151). Middle: a model system with a rigid array of alternating solvent (B) - membrane (A) - solvent 


(B) regions. Right: The Hamaker coefficient, 'H{a,d) in [zj], as a function of the separation between membranes, d, and the 
thickness of the lipid bilayers, a. The Hamaker coefficient is calculated based on the full dispersion spectra of water and lipids 


(hydrocarbons) (Podgornik et al. 20061. 


common additives in paper such as Ti02, and CaCOa in water were obtained as ~ 3 — 19 zJ |Hollertz et al. 2013) 
and can explain important adhesion, swelling and wetting phenomena ubiquitous in paper processing. 


The long-range interaction between proteins is also of vdW nature (Leckband and Israelachvili, 2001 Leckband 


and Sivasankar 1999|) . Estimates for protein-protein interactions across water or dilute salt solutions report Hamaker 


coefficients mostly within the range ~ 10—20 zJ (Farnum and Zukoski 1999). H for interacting proteins, such as 


bovine serum album (BSA), has been found to be ~ 12 zJ by considering a Drude-Lorentz model for the dielectric 
function and the zero Matsubara frequency term included (Neal et al. 1998 Roth et al. 1996). Using the anisotropic 
coarse-grained model of a protein on the level of amino acid residues, one can calculate effective polarizabilities of 
bovine pancreatic trypsin inhibitor (BPTI), ribonuclease inhibitor, and lysozyme in an aqueous solution (Song 20021. 
These results have to be approached with caution, however. Accurate frequency-dependent polarizabilities are rarely 
available either from theoretical or experimental studies (Nandi et al. 2000), thus one has to rely on plausible but 


probably unrealistic model approximations (Song and Zhao 2004). The same is true for the static dielectric constant 
that shows pronounced variation from the inside to the periphery of the protein (Li et al. 2013). The anisotropic 


optical spectrum of collagen, a fibrous protein, has been calculated by ab initio methods (Poudel et al. 2014) and 


used to estimate the corresponding non-isotropic Hamaker coefficients (Dryden et al. 2015). The vdW interactions 


between collagen fibers show a substantial angle-independent component of the Hamaker coefficient ~ 9.3 zJ. The 
origin of the angular dependence of vdW interactions is in fact twofold: the morphological anisotropy, given by the 
shape, and the material anisotropy, given by the dielectric response tensor. Both contribute to the general angular 
dependence and consequently torques between biological macromolecules (]Hopkins et al. 2015), see below. 


B. Lipid membranes 


In general, the vdW interaction is of fundamental importance for the stability of biological matter (Nel et al. 2009) 


and for membrane arrays in particular (Pasichnyk et al. 2008 Petrache et al ., 20061. The essential component of a 


membrane is the lipid bilayer, a planar layer of finite thickness composed of a hydrocarbon core with hydrophyllic 


boundaries facing the aqueous solution (Tristram-Nagle and Nagle 2004). vdW interactions between lipid membranes 


were in fact the first example of using Lifshitz theory in condensed media (Parsegian and Ninham 1969). Calculated 


non-retarded Hamaker coefficients were found to be in the range 1 — 10 zJ. The importance of the ionic screening of 


the zero frequency Hamaker term in electrolyte solutions for membranes has also been carefully quantified (Ninham 


and Parsegian 

1970b 

Petrache et al. 

2006), see above. As already stated, the very high static dielectric constant of 

water ( 

Parsegian 

2006 

1 leads to an anomalously large contribution to the entropy of vdW interactions, which remains 


unretarded for all separations as it corresponds mostly to the n = 0 Matsubara term. However, taking into account 
electrolyte screening at sufficiently large salt concentrations reverses the anomalous effect of the water dielectric 
constant, so that retardation effects emerge from a combination of electrolyte screening and standard retardation 


screening (Ninham and Parsegian, 1970b). 


Most experiments yielding the strength of the non-retarded Hamaker coefficients are actually based on multilamellar 


interaction geometries that allow for detailed osmotic stress small-angle X-ray scattering (SAXS) studies (Kollmitzer 






















































































































et at. 2015 Tristram-Nagle and Nagle 2004) (see Fig. 17). One also needs to consider the non-pairwise additive 


vdW effects in multilamellar geometries that can be significant (Narayanaswamy and Zheng, 2013). The interaction 
of a pair of two lipid membranes with thickness a at a separation d in a multilamellar stack yields for the interaction 
surface free energy density F{a,d): 


J^{a, d) ~ — 


keT 




drr (a + d)^ 

+ ^ f-2^(2 + y, ,) — 2Z(2 + y,l) + Z{2 + y, —^ 
V a + d a + a / 

r> = 1 ' ^ 


(24) 


where ({m,n) is the zeta-function and y = 2^ \/eB(iw„)(a 
with y according to (Podgornik et at. 2006) and A(a;) = ( 


the lipid and water, respectively. Also, p\ b ~ ~ 


(for details see Ref. (Podgornik et al.\ 2006)). 


d). The function Z{2 + y,x) is exponentially screened 
where eA(w) and €5(0;) are the permitivities 
where Q is the magnitude of the transverse wave vector. 


PA^B—PB^A 

\PA^B-\-pB^A 


The non-additive effects in J-{a,d) vanish at d <C a, where the interaction is obviously reduced to that of two 
semiinfinite lipid regions across water. Approximating the water response function by one Debye and twelve Lorentz 
oscillators (Dagastine et al. 2000[ Roth and Lenhoff, 1996| , and the lipid response function by four Lorentz oscillators 
in the ultraviolet regime (Parsegian, 2006) yields H(a = 4 nm, d ~ a) = 4.3 zJ. The usually quoted theoretical 
result with no retardation effects ( jParsegian 2006) is 'H{a^d ~ a) = 3.6 zJ, while experimentally determined H is 
typically in the range 2.87—9.19 zJ for dimyristoyl phosphatidylcholine (DMPC) and dipalmitoyl phosphatidylcholine 
(DPPC) lipid multilayers (Petrache et al ., [1998 ). The lipid bilayer thickness dependence is clearly seen in recent 
experiments with dioleoyl phosphocholine/distearoyl-phosphocholine/cholestrol (DOPC/DSPC/Chol) mixtures that 
phase separate into liquid-ordered (Lo) and liquid-disordered (Ld) domains with H = 4.08 zJ for Ld and H = 4.15 zJ 


for Lo domains (Kollmitzer et al.\ 2015). 


Other fluctuation-induced Casimir-like interactions (Kardar and Golestanian 1999) are also of relevance in the 
context of lipid membranes. Among these, the Helfrich interactions due to steric repulsion between fluctuating 


membranes have received very detailed attention, see (Freund 2013 Lu and Podgornik, 2015). However, more directly 


related to the Casimir effect are the thermal height fluctuations of membranes, constrained on average to be planar, 
that can couple in various manners to the local membrane composition. For example, embedded macromolecules. 


such as proteins (Phillips et al. 2008), cause modifications of the height fluctuations due to the spatial variation of 


in terms of the membrane height function /i(x) (Lipowsky 


the effective membrane rigidity at the position of the inclusion. Thus, there is an elastic Hamiltonian (Deserno, 2015) 

[TM| ) 


H[h(x)] = / 


^Kr{x) (V^/l(x))^ -I- Kr{x) 


d^hix) d^h{y 
dx'^ 


d^h{x)' 

dxdy 


(25) 


where the local bending rigidity Kr(x) and the local Gaussian rigidity Kr(x) are position dependent upon x = 
(x, y) denoting the in-plane coordinates for the membrane projected area (Dean et al.. 2015). In single component 
membranes, where Kr and are constant, the Kr term is zero when the membrane is a free-floating sheet, since by 
virtue of the Gauss-Bonnet theorem it only depends on the boundary and topology of the membrane (David 2004| ). 
The height correlator can be found analytically in this case (Deserno 2015). In the presence of elastic inclusions, both 
moduli contain an unperturbed constant part, Ko,r and Ko,r, as well as the position dependent parts AKr(x), AKr{x') 
that vanish outside of the inclusions. Other parametrizations of the effect of inclusions are also possible and have been 


considered (Netz 1997). Alternatively, membrane inclusions can be also considered as curvature sources (Dommersnes 


and Fournier 1999a). To the above energy one can also add surface energy and an external potential energy when 


appropriate (Noruzifar et al. 2013). Apart from the order of derivatives in the fluctuating field, second for membranes 
and first for electrostatic field, and the local bending rigidity taking the role of the local dielectric permittivity, H[h{x)] 
is completely analogous to the electrostatic field Hamiltonian (Naji et al. 2013) or indeed to the Hamiltonian of critical 


mixtures (Trondle et al.. 2010). Thus one can expect that thermal fluctuations effects will also be present. 

Indeed, inclusions in the membrane, which modify its local mechanical properties, can experience fluctuation- 


induced forces between them I 

Bartolo and Fournier 


Golestanian et al. 1996a|b 

Goulian et al. 1993 Lin et al. 

2011 

Park and Lubensky 

1996 

Yolcu et al. 

2011) 

. The interactions are energetically of the order of ksT, and can 


in certain circumstances, in particular for tensionless membranes, be long-ranged and therefore potentially exhibit an 


important effect on the organization of the membrane (Machta et al. 2012). Several studies have considered coupling 



















































































































































38 


of membrane inclusions to the membrane curvature via the elastic stress (Bitbol et al. 


et al.^ 2011) or to topological defects in orientational disorder (Golestanian et al. 1996a 


2010 Lin et al.. 2011 Yolcu 


Korolev and Nelson, 2008). 


Assuming that Kr(x) and Kr(x) are small, one can calculate the cumulant expansion of the partition function (Goulian 


et al. 1993). The effective two-body interaction between regions deviating from the background rigidity /to,r is then 


obtained in the simple form 


H. = 


keT 




,2 , AKr{x.)AKr{x') 


^'14 


(26) 


after expanding to the lowest order in the deviation from the constant values of the rigidities, A7tr(r) and AKr{r'), 
that vanish outside of the inclusion. When the separation between local regions (inclusions) characterized by change 
in rigidities is much larger than the size of the regions, the first order term in a multipole expansion of the energy 
between the two regions therefore decays as the fourth power of separation. Notably, both ^^(x) and ^^(x), have 
to be present in order to have a fluctuation interaction. On the other hand, considering membrane inclusions as 
curvature sources one can bypass these constraints, at the same time also strongly enhancing and increasing the range 


of the interactions (Dommersnes and Fournier 1999aI. These Gasimir-like forces are dominated by fluctuations and 
their variance also shows a characteristic dependence on the separation ( ]Bitbol et al. 2010) as well as pronounced 
many-body aspects ([Dommersnes and Fournier 1999b), as expected for Gasimir-like interactions. 


Scattering methods developed for the electromagnetic Gasimir effect (Rahi et al. 2011) and discussed in Sec. V 


have been employed for the interaction between two membrane embedded disks discs of radius R to all orders, leading 
to asymptotic result for large separations (Lin et al. 2011). The effective field theory formalism also affords an 
efficient framework for the computation of membrane fluctuation-mediated interactions ( Yolcu and Deserno[ 2012| 
Yolcu et al.. 2011). In the case of broken cylindrical symmetry of the inclusions, the fluctuation interaction retains the 
same separation dependence, but its strength depends on the two orientation angles as cos 26i cos 202 and cos 4(0i +02) 


Golestanian et al. 

1996b 

Park and Lubensky 

1996 


This of course implies the existence of fluctuation or vdW 

torques, see below. 

Experimentally, fluctuation mediated interactions between membrane inclusions might be difficult to measure di¬ 
rectly, if at all feasible. More promising seems to be the detection of their consequences, like fluctuation-induced 


aggregation of rigid membrane inclusions (Dommersnes and Fournier, 1999b 


the miscibility of lipid mixtures in multicomponent membranes (Dean et al. 


Weikel, 2001) or through their effect on 


2015 Machta et oL 2012). 


C. van der Waals torques 


Biological materials are typically anisotropic in terms of shapes as well as response properties (Hopkins et al. 


2015). Such anisotropy leads to vdW torques, which have been first studied by Kats (Kats 1978) for the special case 


of isotropic boundaries with anisotropic intervening material, and independently by Parsegian and Weiss (Parsegian 


and Weiss 1972) who studied the inverse case of bodies with anisotropic dielectric response interacting across an 


isotropic medium in the non-retarded limit (Kornilovitch 


much later in a veritable tour de force by Barash ( Barashj 


2013). The full retarded Lifshitz result was obtained only 


1978), following previous partial attempts (Barash 1973), 


leading to a series of recent developments (van Enk 19951 Munday et al.. 2005, 2008 Shao et al. 2005). The general 


Lifshitz formulae for the interaction between two anisotropic half spaces or even an array of finite size slabs (Veble and 


Podgornik 2009) are algebraically very complicated and untransparent, with little hope of a fundamental simplification 


et al. 

2009 

2013c 

2014) 


(Philbin and Leonhardt 2008). Morphological anisotropy effects are seen either between anisotropic bodies (Emig 


Rahi et al.. 2011) or between surfaces with anisotropic decorations such as corrugations (Banishev et al. 


2013c 2014). vdW-like torques have been predicted also between anisotropic inclusions within fluctuating membranes 


(Golestanian et al. 1996b Park and Lubensky 1996). 


The first attempt to evaluate the vdW interaction between two cylinders comes from Barash and Kyasov (Barash and 


Kyasov 1989). Results for two infinitely long anisotropic cylinders can be obtained in a dilution process (Parsegian 


2006), such that the presence of dielectric cylinders can be considered as a small change of the dielectric permittivity 


of two semi-infinite regions (Pitaevskii 2008). This approach leads to interactions between infinite cylinders of radius 


R at minimum separation d, at any angle of inclination 0 in non-retarded (Rajter et al. 2007) as well as retarded 


limits (Siber et al. 2009), but also between cylinders and anisotropic semiinfinite layers (Hopkins et al. 2015 Saville 


et al. 2006). The vdW interaction free energy for inclined cylinders is 

(^i?2)2 


Rid, 9) = - 


27r d'^ sin 0 


(M(°)((i)+M(2)(‘^)cos20) , 


(27) 
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HdJ) = (^('»(d) + ^P)(d)(2-cos25)) 



FIG. 18 (Color Online) ) Forces and torques between anisotropic molecules/molecular aggregates. Top: f{d, 6), torque between 
two inclined long cylindrical molecules (above DNA below SWCNT) as a function of the closest separation d and angle of 
inclination of their cylindrical axes 9, calculated within the Lifshitz theory. The two Hamaker coefficients A’'°^((i) and A^^\d) 
depend on the dielectric anisotropy of the interacting materials as well as on the separation due to retardation effects, but they 
do not depend on mutual orientation. Bottom: r(d, 0), torque per unit surface area between two planar semi-infinite slabs at 
separation d with angle of inclination 9 between their anisotropic uniaxial tensors. Only the surface layer of the bottom slab 
and one molecule of the upper slab are shown. Again, the Hamaker coefficient A^^\d) depends on the dielectric anisotropy 
of the interacting materials as well as on the separation due to retardation effects, but not on the orientation of the dielectric 
tensors. Adapted from (Dryden et al. 20151. (Graphics courtesy J.C. Hopkins). 


wher e 1/ sin 6 stems from the shape anisotropy and the cos 26 dependence associated with the material anisotropy. 
Fig. 


18 


The Hamaker coefficients, and are functions of separation and the relative anisotropy 

measures, but do not depend explicitly on the angle of inclination 9. They also depend on the material types and 


Matsubara sampling frequencies (Siber et am200^ [Stark et al. 20151. In the symmetric interaction case (Hopkins 


et al. 20151 both ^(°^((i) and A^'‘^{d) decompose into a square and thus cannot be negative or change sign. In the 
asymmetric case, however, the sign of the interaction as well as the sign of the torques are more complicated, as 
they depend on the perpendicular and parallel dielectric response of the interacting bodies. They do not follow the 
general rule for interacting planar bodies 1 and 2 across a medium m, with a change of sign implied by the sequence, 
£1(0;) > £m(w) > S 2 {uj), see above. 

vdW torques between semi-infinite anisotropic materials also imply torques between anisotropic cylindrical 


molecules, such as filamentous graphitic systems of metallic and semiconducting single-walled CNTs ( Rajter et al. 
2013, 2007[ Siber et al.. 20091, multiple composites of DNA (Sontz et al. 2012 Young et al. 2014[ Zheng et al.. 20091, 


type I collagen (Cheng et al.. 20081, and polystyrene (Jin et al. 20051. DNA optical properties (^Pinchuk 2004) have 


been converted into the corresponding separation dependence vdW free energy in the case of pairs of nucleotides 
(Pinchuk and Vysotskii 1999 2001). A static dielectric constant of 8 for DNA, that enters the n = 0 Matsubara 
term, was recently measured inside single T7 bacteriophage particles by electrostatic force spectroscopy (Cuervo 


et al. 2014). Optical dispersion data are also available for single nucleotides, nucleosides and derivatives, synthetic 
polynucleotides (polyuridylic acid, polyadenylic acid, poly-AU), various nucleic acids, such as RNA and native bacte¬ 
rial DNAs in aqueous solutions (Voet et al. 1963), or wet and dry polymerized oligonucleotides and mononucleotides 


(Silaghi et al. 2005 Zalar et al. [2007 ) as well as dry DNA thin films (Sonmezoglu and Sonmezoglu 2011). Optical 
properties of DNA oligonucleotides (AT)IO, (AT)5(GC)5, and (AT-GC)5 using ab initio methods and UV-Vis decadic 
molar absorbance measurements show a strong dependence of the position and intensity of UV absorbance features 
on oligonucleotide composition and stacking sequence ( [Schimelman et oT 2015). The calculated Hamaker coefficients 
for various types of DNA molecules are overall small but depend on the base-pair sequence details and could control 


the finer details of the equilibrium assembly structure (Bishop et al. 2009). In fact, the stacking sequence dependence 


of the optical properties has important repercussions for the molecular recognition between two approaching DNA 
molecules that depends on vdW interactions (Lu et al.. 2015). The angle-independent part of the Hamaker coefficient 
is ~ 5zJ, while the angular part is effectively zero when the zero-frequency Matsubara component is fully screened 
by the electrolyte solution. At least for DNA molecules, it then appears that the anisotropy effects stem purely from 
the shape anisotropy, whereas this is not the case for CNTs. Among the fibrous proteins collagen also shows strong 
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vdW interactions with silica resulting in a Hamaker coefficient that is 39% larger than that of the silica-(GC)10 DNA 
interaction at 5 nm separation ( Dryden et aH['2015 ). 

Though the vdW torque, define d as T((i, 9) = - eminently measurable, it has however not been measured 

directly yet (Capasso et al. 2007 Chen and Spence 2011 lannuzzi et al. 2005). The anisotropy that engenders the 
vdW torque can be of different origins: it can result either from anisotropy of the dielectric response of the interacting 


bodies or from their asymmetric shape (Hopkins et al. 2015). Both give effective Hamaker coefficients that depend 


on the mutual orientation of the dielectric or shape axes (Dryden et al.. 2015). The material anisotropy can also be 
either intrinsic, or a consequence of arrays of nanoparticles embedded in an isotropic background (Esquivel-Sirvent 


and Schatz 2013), and/or a consequence of the action of external fields (Esquivel-Sirvent et al. 2010). It remains 


unclear which anisotropic effects would be best suited for accurate experiments. 


D. Electrostatic fluctuations 


The n = 0 (’’static”) Matsubara term corresponds to the classical partition function of the Coulomb system and 
can have a form very different from the one derived from the Lifshitz theory. For an interacting Coulomb fluid 
such as a confined electrolyte or plasma, a counterion only system, or a system of dipoles or polarizable particles in 
an inhomogeneous dielectric background, the n = 0 Matsubara term corresponds to the free energy of fluctuations 
around the mean-field in the range of parameter space where the mean-field (Poisson-Boltzmann or weak coupling) 
approximation holds (for details, see (Naji et al. 2013| )). This leads to effects such as screening of the Hamaker 
coefficient, universal value for the Hamaker coefficient or its anomalous separation dependence. In effect, the n = 0 
Matsubara term actually corresponds to Caussian or one-loop electrostatic potential fluctuations around the mean- 
field for a fully coupled system ( |Netz 2001b Podgornik and Zeks, 1988) and thus presents a first order correction to 
the description of Coulomb fluids on the mean-field level (Holm et al.. 2001). Formally this follows from an effective 


non-Caussian “field-action” that one can derive from an exact field-theoretic representation of the confined 


Coulomb fluid partition function in terms of the fluctuating local electrostatic potential (Edwards and Lenard 1962 


Podgornik 1989). The mean-field theory is then defined as the saddle-point of this field action (Naji et al. 2013). 


Thermal fluctuations around the saddle-point at the first-order loop expansion representing the contribution from 


correlated Caussian fluctuations around the mean-field or saddle-point solution (Dean et al. 2014), lead to a thermal 


fluctuation-induced vdW-like attraction in the form of the trace-log of the “field-action” Hessian 

S^S[^] 


F = —ksT TrLog 


i5^(r)d(/(r') 


0 (</"), 


(28) 


0(r)=-iV>pB(r) 


where ^(r) = —i'!/pB(r) is the saddle-point (Poisson-Boltzmann) electrostatic potential configuration. This second- 
order correction universally lowers the interaction pressure between surfaces and thus leads to an attractive contribu¬ 
tion to the total interaction pressure. It also includes a term that exactly cancels the zero-frequency contribution in 


the Lifshitz theory (Podgornik 1989) so it should be viewed as a substitute for the n = 0 Lifshitz term. As a rule, this 


fluctuation attraction is weaker than the repulsive leading-order saddle-point contribution, thus the total interaction 


remains repulsive (. 

Metz 2001a 

Netz and Orland 

20001. In certain models that either assume charge asymmetry 

(Kanduc et al. 

o 

o 

1 ), surface condensation or adsorption of counterions on (hxed) charged boundaries, the repulsive 


mean-held effects are strongly suppressed and the total interaction is then mostly due to thermal huctuations (Lau 


et al. 

2001 

2002 

Lau and Pincus 

2002 


This can happen in the case of oppositely charged surfaces where the 


1999D . 

A related problem of thermal electrostatic huctuations is presented by the Kirkwood-Shumaker (KS) interactions 


that exist between macroions with dissociable charges, such as proteins ( 

Lund and Jonsson 

2013). Originally this 

interaction was obtained from a perturbation theory around an uncharged state ( 

Kirkwood and Shumaker, 1952 

Lund 


and Jonsson 2013). The KS interaction is similar to the thermal vdW interaction but it corresponds to monopolar 


charge huctuations (Adzic and Podgornik, 2014) and is thus in principle much longer ranged. Monopolar huctuation 


cannot arise for hxed charges on interacting bodies and some surface charging mechanism or charge regulation, where 
the macroion surfaces respond to the local electrostatic potential with a variable ehective charge, is needed (jBorkove'c 


et al.\ 2001). Formally charge regulation can be described by another, nonlinear source term in the held-action. 


/ 5 ((/(r)), that involves the huctuating potential at the surface (S) (Adzic and Podgornik 2015). Nonlinearity of this 


held-action is essential as a linear dependence on the huctuating potential, in fact, corresponds to a hxed charge 
that cannot exhibit monopolar charge huctuations. These are given by the surface capacitance determined by the 
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second derivative of /s(0) with respect to the surface potential, and the KS interactions depend quadratically on 
this capacitance. The KS interaction between two macroions in the asymptotic regime between particles 1 and 2 
then assumes the form T ~ The exact form of /s and thus the capacitance C is not universal as it depends 


on the surface-ion interaction (Adzic and Podgornik 2015 Fleck and Netz, 2007 Markovich et al.. 20141. Although 


calculating the KS interaction is challenging, exact solutions are available in ID, demonstrating a rich variety of 
behaviors due to charge regulation and the ensuing correlated fluctuations (Maggs and Podgornik 20141. Some of 


these ID properties transfer also to the more realistic 3D models between globular proteins with dissociable surface 
charge groups ( [Adzic and Podgornik 2015). Similar anomalously long-range monopolar fluctuations and concurrent 
Casimir/vdW-like interactions can also result from a different mechanism where monopolar charge fluctuations result 
not from charge regulation but rather from nano-circuits with capacitor components, where fluctuating charges are 


transferred through the wire connection in a capacitor system (Drosdoff et al. 2015). 


The link between Coulomb interactions and thermal Casimir/vdW interactions has been implicated also in some 
theoretical approaches to the Hofmeister or specific ion effects (Salis and Ninham 2014). These works motivated 


numerous investigations of non-electrostatic ion-specific interactions between ions and surfaces and their role in mod¬ 
ifying surface tension of electrolyte solutions or indeed the solution behavior of proteins. In fact, the standard 
Onsager-Samaras result was only recently realized to be fluctuational in nature (Markovich et al. 2014). Ninham and 
coworkers as well as others (Edwards and Williams 2004) made attempts to include vdW interactions into a complete 


theory of ion interactions in confined aqueous solutions (Bostrom et al. 2006 Bostrom and Ninham 2006 Bostrom 


et al.. 2005). The major problem in including vdW interactions into the description of inhomogeneous electrolytes is 
that the contributions of fixed charges and ion polarizability are in general not additive ( [Demery et al. 2012). How¬ 
ever, they sometimes can be approximated by an additive contribution on the strong coupling level, provided that 
the polarizability of the ions is large enough. A popular Ansatz that simply adds a vdW ion-polarizability dependent 
contribution to the electrostatic potential of mean force has thus a very limited range of validity. 


VII. EXPERIMENTS PROBING MATERIALS ASPECTS OF VDW/CASIMIR INTERACTIONS 


The theory and computation of vdW/Casimir and related fluctuation-induced interactions are experiencing an ex¬ 
pansion paralleled by numerous recent materials and structured systems discoveries. Experimental efforts have also 
been reported probing how different materials can be utilized to modulate this subtle force. Most recent experiments 
have concentrated primarily on structured systems and some biological matter, as discussed in this review. Addition¬ 
ally, novel measurements have served as a propeller to the field not only to validate certain theoretical predictions, 
but also to identify new problems. 

Atomic force microscopy (AFM) techniques have been employed in vdW measurements giving unprecendented 
insights into smaller and hetergoeneous systems. Pulling single molecules with an AFM tip from metallic surfaces 
(Wagner et al.. 2014) confirms the asymptotic d~^ force law and quantifies the non-additive part of the vdW interac¬ 


tion. Non-additive effects have also been demonstrated in adhesion measurements in various tribological environments. 


as well (Loskill et al. 2012 2013). Another recent report gives a clear evidence of the vdW screening capabilities in 


graphene/MoS 2 heterostructures giving insight into adhesion properties of graphene and other layered materials (Tsoi 
et aL\ 2014). These non-additive and screening effects are particularly challenging for theory, which has motivated 


developments in first-principles calculations methods, as discussed in Sec. II. 


The first experimental measurement of the Casimir force involving graphene has been reported in (Banishev et al. 


2013a). Good agreement with theory taking into account the Dirac spectrum has been achieved for the studied 
graphene/Si02 setting (Banishev et al.. 2013a Klimchitskaya and Mostepanenko 2015). However, many issues 


need to be investigated further. For example, more precise comparison with theory is needed, which on the other 
hand requires measurements of free standing graphene interactions. Schemes to determine the asymptotic distance 
dependence and temperature effects are absent. Experiments in this direction will be extremely desirable as they can 
serve as a validator for the numerous theoretical predictions. Reliable experiments for stacks of graphene are also 
much desired in order to determine directly the binding energy of graphite and settle the wide range of values reported 


through indirect measurements (Benedict et al. 1998 Zacharia et al. 2004). It would also be very interesting to seek 


experimental knowledge of vdW/Casimir interactions involving other systems with a Dirac spectrum. Probing novel 
topological phases in such dispersive interactions using topological and Chern insulators can be very beneficial. 
Experimental measurements in structured materials continue to shed light on many aspects of the Casimir force. 


from non-additive effects in grating and related geometries (Banishev et al. 2013c Chan et al. 


2013), to repulsive interactions in interleaved structures (Rodriguez et al. 2008a Tang et a/.[|2015[ ), to strong tem¬ 
perature corrections arising at large separations (Sushkov et al. 2011) or in situations involving structured magnetic 
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media (Bimonte 20151. Experiments at nano-metric scales, involving objects with non-trivial surface topology due 


to roughness or patch charges, are also beginning to push the boundaries of theoretical techniques (Kim et al. 2010 


Lamoreaux 20101. For instance, recent AFM characterizations of the surface morphology of planar objects (van Zwol 


et al.\ 20081 coupled with predictions based on the above-mentioned state-of-the-art simulation techniques (Broer 


et al.i 2012 j reveal that the presence of roughness on the scale of their separations manifests as strong deviations in 


the power-law scaling of the force (Broer et al. , 2013). Challenges that continue to be addressed in current-generation 


experiments include the need to control and calibrate materials properties. Specifically, accurate comparison between 
theory and experiments require accurate knowledge of the dielectric response over a broad spectral region, spurring 


recent efforts to characterize numerous material properties from DC to ultraviolet wavelengths (Sedighi et al. 2014 


van Zwol et al. 2009). Finally, a number of experiments employing novel fabrication techniques have begun to explore 
Casimir forces in on-chip, integrated systems (Yamarty and McNamara 2009 Zho et al.. 2013|) where parallelism is 


no longer a key impediment, removing the need for external instrumentation needed to bring objects closer together 
and paving the way for applications where the force is exploited in conjunction with other effects (e.g. mechanical or 


2009). 


optical actuation) to enable new functionalities (Pernice et al. 2010 Rodriguez et al.. 2011b Yamarty and McNamara 


On the other hand, small-angle X-ray scattering techniques are useful for measuring vdW potentials in soft- and 
bio-matter systems. Their contribution to the elucidation of the details of long-range interactions in multilamellar 
membrane array context is crucial and continuing (Kollmitzer et al. 2015). It remains to be seen whether SAXS 


coupled to osmotic stress could provide also a vdW component to the interaction potential between filamentous 


bio-molecules (Yasar et al. 2015), allowing also a determination of the Hamaker coefficient that could be compared 


with calculations (Poudel et al. 2014 Schimelman et al. 2015). While the Hamaker coefficients for general proteins 


could be approaching solid predictions (Eifler et al. 2014), their experimental determination is marred by the limi¬ 


tations inherent in the second virial coefficient determination of the global characteristics of the interaction potential 


(Prausnitz, 2015). Novel methodologies such as colloid-probe AFM could be a new potential source of valuable data 
on intermolecular potential (Borkovec et al. 2012), including protein-protein interactions at all values of separation 


(Singh et al. 2015). 


VIII. FUTURE OUTLOOK 


The Casimir/vdW force has manifestations in many parts of physics as discussed at length in this review. The 
quest for a fundamental understanding of this ubiquitous and subtle force bridges concepts from condensed matter 
and high energy physics, which has become much more apparent with recent discoveries of novel materials. This field 
has also stimulated the development of computational methods at the atomistic levels as well as larger scale with the 
goal of taking into account the collective and non-additive nature of the same dispersive interaction. This particular 
direction has also been stimulated due to materials science expansion and better design of devices. Nevertheless, the 
field can become even broader with several eminent problems awaiting solutions. 

Materials with Dirac spectrum hold much promise to discover new science about the vdW/Casimir interaction. 
Unusual behavior of the Casimir force in terms of sign, magnitude, distance dependence, and other factors has 
been found in graphene, TIs and CIs, as discussed in this review. But there are many new entering players with 
emergent properties. 2D TIs, such as HgTe/Cd/Te, Bi bi-layers, and InAs/GaSb; 3D TIs, such as Bii_xSbx, Bi 2 Te 3 , 
Heusler alloys, and topological crystalline insulators, such as SnTe, Pbi_xSnxSe, have prominent spin-orbit interaction, 
which coupled with the Dirac spectrum can lead to diverse optical response (Welding et al. 2014). 3D Weyl and 


Dirac semimetals, such as Cd 3 Al 2 , Na 3 Bi can also be put in this category. Recent reports show that the quantum 
electrodynamics in Weyl semimetals results in a non-trivial response due to the associated axion field, which can lead 


to a repulsive Casimir interaction (Grushin 2012 Wilson et a/. 2015). However, external electric and magnetic fields 


together with temperature and doping can modulate the electronic structure and optical properties by creating new 
topological phases, such as valley polarized materials for example. These are yet to be studied in the context of the 
vdW/Casimir interaction. 

Another interesting direction to explore originates from nonlocality, especially at very small separations. For 
instance, at nano-metric scales, the combination of structured materials, thermal as well as dielectric inhomogeneities, 
and non-local effects associated with atomic-scale physics, can potentially conspire to affect fluctuation phenomena. 


Preliminary works studying non-local material effects have begun to shed light on these issues (Despoja et al. 2011 


Esquivel-Sirvent et al. 2006 Luo et al. 2014aI. One promising set of scattering techniques that could be used to 


tackle these emerging regimes are volume-integral equation methods, related to surface-integral equations but involving 


volume rather than surface currents inside the bodies (Polimeridis et al. 2014). Another set of techniques that are 
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beginning to pave the way for fundamentally new designs in nanophotonics but which have yet to be exploited in 
Casimir computations are large-scale optimization methods (Bendsoe and Sigmund, 20031. While such brute-force 


explorations require careful and efficient formulations due to the large number of required calculations, the above- 
mentioned numerical developments offer hope that such an approach to design is within reach. Finally, although many 
of the interesting, non-additive Casimir effects predicted thus far still remain out of reach of current experiments, 
perhaps related physical principles can be employed to discover other structures where non-monotonicity and/or 
repulsion is larger and more experimentally accessible. 


Transformation optics, a powerful method for solving Maxwell’s equations in curvilinear coordinates (Leonhardt 


2006 Pendry et al. 2006), may offer a different perspective to the vdW/Casimir effect, especially for systems that 


have sizes comparable to their separation. It can be an efficient numerical approach for vdW calculations by taking 
into account nonlocal effects for absorption and scattering spectra, electromagnetic modes and field enhancement. 
Transformation optics has been a powerful tool for optics design with applications, such as perfect lensing and 
cloaking (Pendry et al. 20061. Recent reports have shown interesting physical insight for vdW interactions in 3D 


objects with nonlocal dielectric properties (Luo et al ., 2014b Pendry et al.\ 2013). An exciting future direction can be 


to examine many of such predictions and applications in the context of vdW/Casimir interactions for new directions 
of control and manipulations. 

We further note that fluctuation-induced phenomena go beyond the dipolar fluctuations that give rise to the 
vdW/Casimir force. Charge and potential fluctuations, beyond the situations discussed in bio-materials, may be very 
interesting in solid state devices. Dispersive forces of charged objects are much less studied. Systems with reduced 
dimensionality may be used to investigate much longer ranged monopolar fluctuation forces which can exist on their 


own or be entangled with the ’’traditional” dipolar fluctuations (Bimonte 2007 Drosdoff et al. 2015). This practically 


unexplored direction holds promise to expand fluctuation induced interactions beyond dipolar excitations and further 
broaden the perspective of Casimir-like phenomena. 

We also would like to mention the puzzle about the relaxation properties of conduction carriers and their role in the 
Lifshitz theory. Depending on the dielectric model used, Drude-like or plasma-like, different magnitudes of the thermal 
Casimir force between metallic or magnetic systems are predicted. The crux of the issue is the correct description of the 


low-frequency optical response of the materials. In (Sushkov et al.\ 2011) the Casimir force between metallic samples 


was measured, and the authors interpreted their results in agreement with the Drude model, after subtracting a force 
systematics due to electrostatic patches that was modeled and fitted to the total observed force. Recent independent 
measurements of patch distributions on metallic samples used in Casimir force experiments report different strengths 
and scaling laws for the patch contribution to the total force ([Behunin et al. 2014 Garrett et al. 2015). In contrast, 


several other Casimir experiments (Banishev et al. 2013b Decca et al. 2005) seem to be in agreement with the 
plasma model description, which is surprising given that this model neglects dissipation effects in metals. More 
recently, a proposal was put forward ( ]Bimont^ 2014a|b ) based on the isoelectronic technique (which eliminates the 
need for electrostatic corrections due to patches) that makes a significant step forward towards a resolution of this 
controversial issue. With this set up, it becomes possible to strongly enhance the discrepancy between the predictions 


for the Casimir force based on either model for the dielectric response. Preliminary measurements (Decca, 2015) are 


in favor of theoretical extrapolations to low frequency based on the plasma model. To date, a basic understanding 
of this fundamental problem in Casimir physics is still missing. Perhaps one way to resolve this issue is to consider 
materials with reduced dimensions and novel phases. For example, the thermal fluctuations effects are much more 
prominent for graphene, which can be a possible direction to explore in this context. A different pathway could be 
that more sophisticated models for the response properties are needed. Nevertheless, a possible resolution to this 
open problem may be found by improving our understanding of materials properties. 


IX. CONCLUSIONS 

A broad perspective in materials and their properties was given to the field of van der Waals and Casimir inter¬ 
actions. This comprehensive review shows that this a broad area where materials have played an important role in 
motivating the development of new theoretical models and computational approaches as well as advances in exper¬ 
imental techniques. Materials may hold the answers of several open problems in fluctuations-induced phenomena, 
which are of fundamental and applications relevance. 
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